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Abstract of the Dissertation 
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Measurement of Josephson- Junction 

Qubits 

by 
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Doctor of Philosophy 
in 
Physics 
Stony Brook University 
2005 



This work deals with two pressing issues in the design and opera- 
tion of Josephson qubits - loss of coherence and measurement. 

The first part is devoted to the problem of decoherence in single 
and double qubit systems. We present a quantitative description of 
quantum coherent oscillations in the system of two coupled qubits 
in the presence of weak noise that in general can be correlated 
between the two qubits. It is shown that in the experimentally 
realized scheme the waveform of the excited oscillations is not very 
sensitive to the magnitude of decoherence correlations. Modifica- 
tion of this scheme into a potentially useful probe of the degree of 
decoherence correlations at the two qubits is suggested. Further- 
more, an explicit non-perturbative expression for decoherence of 
quantum oscillations in a single and double qubit systems by low- 
frequency noise is presented. Decoherence strength is controlled by 



iii 



the noise spectral density at zero frequency while the noise correla- 
tion time r determines the time t of crossover from the l/y/i to the 
exponential suppression of coherence. The first part is concluded 
by introduction of a useful Monte-Carlo-based numerical technique 
for simulations of realistic qubit dynamics subject to general noise, 
time dependent bias and tunneling, and general environment tem- 
perature. 

The second part addresses the problem of quantum measurement 
of qubit through discussion of a versatile new type of the magnetic 
flux detector which can be optimized to perform several measure- 
ment schemes in quantum-limited regime. The detector is based on 
manipulation of ballistic motion of individual fluxons in a Joseph- 
son transmission line (JTL), with the output information contained 
in either probabilities of fluxon transmission/reflection, or time de- 
lay associated with the fluxon propagation along the JTL. The 
equations for conditional evolution of the measured system both 
in the transmission/reflection and time-delay regimes are derived 
for this detector. Combination of the quantum-limited detection 
with control over individual fluxons should make the JTL detector 
suitable for implementation of non-trivial quantum measurement 
strategies, including conditional measurements and feedback con- 
trol schemes. 
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Chapter 1 
Introduction 



Any algorithmic process can be simulated efficiently (i.e. in polynomial time) 
using a Turing Machine is the celebrated Church- Turing thesis that for a long 
time instilled the belief that the efficiency of computation is independent of its 
underlying physical process P1I2]- F° r this reason, the efficiency of computer 
architecture built on irreversible binary language was not questioned at first. 
Rather, one of the important topics in the earlier research of computational 
devices was the energy dissipation. One such study [3] showed that the erasure 
of information is the only source of energy dissipation. Later, it was suggested 
[I] that the fundamental limit of energy dissipation, ksTln2 per one erased 
bit of information, can be avoided by employing the computation that is both 
physically and logically reversible since reversible system can always be "rolled 
back" to the zero state along constant energy path. 

A particular way to employ reversible computation is by using quantum 
mechanics as the underlying physical process of computation since all closed 
quantum systems are reversible. In addition to the vanishing energy dissipa- 
tion, a quantum system provides a possibility of preparing a superposition of 
its different states. This speeds up the information processing since it allows 
operating on many states of the computational basis simultaneously. 

The advantage of this so-called quantum parallelism and the reversibility 
of quantum systems was for the first time fully realized by David Deutsch. 



2 



He formulated quantum computation as a new field of computation based on 
quantum rather than classical dynamics [S]. The same work featured the first 
quantum algorithm that was in certain cases more efficient than any Turing 
machine algorithm for the same problem. In 1992 together with Jozsa, Deutsch 
improved on this algorithm so that it was always more efficient than any clas- 
sical one jH] . Further discoveries of quantum algorithms that are more efficient 
than the respective classical algorithms by Simon jjj, Shor [H], and Grover jU] 
propelled quantum computation into the forefront of scientific research. 

The fundamental element of quantum computation is the quantum bit 
or qubit for short. It is a two-level unitary system where the unitarity is 
manifested through the presence of the quantum coherent oscillations. The 
oscillations represent periodically evolving probability of a qubit occupying 
some mixture of its two states. For this reason reconstructing previous states of 
the qubit is always possible. Decoherence, i.e. the loss of coherence, suppresses 
the quantum coherent oscillations and with it the reversibility of the qubit and 
in the process degrades its ability to preform the tasks needed in quantum 
computation [TU] . 

An isolated quantum system does not exhibit decoherence. For it to hap- 
pen, the quantum system needs to interact with other systems in a such way 
that its dynamics becomes dependent on the surrounding systems. Classi- 
cal systems in such circumstances undergo energy relaxation as described by 
the fluctuation-dissipation theorem [TT]. Quantum systems exhibit energy re- 
laxation and decoherence. While quantum relaxation has its counterpart in 
classical dynamics, decoherence is a purely quantum phenomenon. It is caused 
by quantum phase interference or entanglement of two systems. It is not con- 
ditional on the systems being in contact all the time, rather if entangled in 
the past, the two systems generally stay entangled. 

For computational purposes the qubit needs to be initialized, coupled with 
other qubits, and in the end measured. This requires entanglement of the 
qubit with the other qubits and the surroundings. Thus to have a qubit 
completely isolated from the environment is impossible. Rather, it is necessary 
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to find ways to design physical qubits so that they exhibit prolonged coherent 
oscillations while integrated in larger structures. 

The leading candidate for achieving such goals are superconducting Joseph- 
son junctions. They are scalable and integrable solid state devices where 
Cooper pairs of electrons move in phase and behave as a single quantum 
particle of much larger size than than the elementary quantum particles. Ad- 
ditionally Josephson junctions are ideal for nano-fabrication and integration 
into larger circuits needed for additional tasks such as measurement. 

Measurement can be seen as an extraction of information from the com- 
putational device and its transcription into a classical signal. The rate of the 
information extraction from the qubit is at the most as fast as its rate of deco- 
herence due to the back-action of the measuring device ^2] • Furthermore, the 
"no-cloning" theorem (eg. see prevents us from making multiple copies 
of the same qubit and repeating the same measurement to recover any lost 
information. Therefore, when devising the measurement of a qubit it is im- 
perative to be as close as possible to this ideal limit, otherwise the non-ideal 
measurement of qubit will result in the loss of information. 

We see that decoherence is present and measurement is needed in the design 
of quantum computers. Understanding these processes is imperative for the 
further advancement of the field. 

1.1 Quantum Mechanics of a Two-Level Sys- 
tem 

We can think of a qubit as a quantum system where two of its neighboring lev- 
els with similar energies are separated by big energy gaps from the remaining 
levels, so that they can be considered independent. This can be accomplished 
with a mesoscopic Josephson junction consisting of a small superconducting 
island or a box separated by a thin insulating barrier from a larger supercon- 
ducting electrode so that Cooper pairs can tunnel to and from the box [T4*l li5|. 
If the superconducting gap of the box is much larger than the thermal energy 



4 



ksT and the charging energy of the box E c = 2e 2 /C, then its dynamical prop- 
erties are solely determined by the excess number of Cooper pairs in the box. 
The potential difference V q between the island and the electrode induces con- 
tinuous polarization charge CV q inside the Cooper pair box of capacitance C. 
Because of the small size of the box, its capacitance is small, and the energy 
spectrum of the box is dominated by its charging energy. The quantum tun- 
neling of Cooper pairs between the box and the reservoir is then a perturbative 
effect described by Josephson tunneling amplitude —Ej/2. 

In the basis spanned by the number of Cooper pairs in the box \n), the 
Hamiltonian of the system consists of the charging and tunneling parts which 
in that order yield the Hamiltonian as: 

H = E c (n - n q f - ^ (\n)(n + l| + \n + l)(n\) , (1.1) 

where n is the excess number of the Cooper pairs in the box, and n q = CV q /2e 
is the induced static charge expressed as the number of Cooper pairs, n q e 
[0, 1). If the applied potential V q is adjusted that n q « 1/2, the gap between 
the system's two lowest charge states \n = 0) = \R) and \n = 1) = \L) is 
small (i.e. E w = Ec(l — 2n q ) m 0), while the gap to the other states is on 
the order of Ec and considerably larger. If k^T <C Ec, the two lowest levels 
are separated from the rest of the levels (figure II. lb ) and as experimentally 
demonstrated in [THIE] they can be considered independent. 

After truncation of the higher energy levels, the Hamiltonian (jl.lj) can be 
written in the form of a general qubit Hamiltonian 

H = --{ea t + Aa x ), (1.2) 

where a x>z are Pauli matrices, e = Ec{l — 2n q )/2, and A = \Ej\ 1 . The 

eigenvectors |0) and |1) of this Hamiltonian define the energy basis of the 

qubit with eigenenergies ±0/2. The qubit Hamiltonian is diagonal in this 

the case of a two-level system, an individual A can always be made real by a unitary 
transformation. 
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Figure 1.1: a) Energy spectrum of Josephson junction operating in charge 
regime with dots representing energy values of different charge states and 
dashed lines representing two lowest quantum levels with charging energy de- 
generacy removed by Josephson tunneling (Eq/Ej — 8, q— 1/2). b). Picture 
and scheme of the first operational charge qubit. c) The quantum coherent 
oscillations observed by Nakamura et. al. with the device pictured above fH] 
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basis, H diag = 



—Q/2 a z , and the vectors spanning the two bases are related as 



\L) = ar_|0> +x+\l), \R) = x+\0) - x-\l) 



(1.3) 



where x± = ^{Vt±e)/2Vt and Vt = Ve 2 + A 2 . 

The state of the qubit is fully specified by its wave function \ty(t)), while 
its time evolution is determined by unitary transformation on 



Given the normalized wave function of the qubit as 1^(0)) = a|0) + 6|1), in 
time t it evolves to |^(t)) = ae tQt ^ 2 \0) + be~ int / 2 \l) expressed in the energy 
basis, or equivalently to 



¥(*)) = {ax + e int / 2 - bx.e- mt ' 2 ) \L) + {ax.e 1 ^ 2 + bx + e~ int / 2 ) \R), (1.5) 



if expressed in the charge basis. As it can be seen above, the probability of 
finding the qubit in one of its eigenstates is stationary, while the probability of 
finding the qubit in one of its charge states (jl.5|) exhibits coherent oscillations 
with period 2ir/Q. 

Quantum coherent oscillations are signatures of both the parallelism and 
reversibility of quantum systems. Thus, the observation of coherent quantum 
oscillations in any qubit design is imperative. The time of coherence expressed 
as the multiple of time interval that takes to perform a single operation on a 
qubit define the quality factor of a qubit. 

The first successful solid state qubit design ^H] was a Cooper pair box 
qubit ( figure ITTTb ). It has demonstrated both the energy band structure of 
Cooper pair box and the existence of the quantum coherent oscillations of 
Cooper pairs ( figure ITTTb ). The quality factor of the qubit was low due to 
the decoherence induced by the proximity of the measurement electrode. The 
subsequent experiments |19|I2U[ I2T] have improved on the shortcomings of the 
early experiments, but the same experiments have labelled the suppression 




(1.4) 
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of decoherence as the most pressing problem in building a qubit suitable for 
realistic quantum computation. 

1.2 Josephson Junctions 

The Cooper pair box introduced above is an example of a Josephson Junction 
that in general refers to any two superconducting electrodes separated by a 
weak link and characterized by the presence of a zero- voltage supercurrent as 
postulated by Brian David Josephson in 1962 [221 122]- The supercurrent I s 
across the junction is related to the difference A0 between the wave function 
phases of Cooper pair condensates of the two electrodes: 

I s = J c sin(A(/>), 

where I c is the critical current. Zero voltage persists across the junction as 
long as the total current across it is smaller than the critical current. In 
this case, the phase difference remains stationary and the junction is in S- 
state. Increasing the current above the critical current switches junction into 
a running or R-state, finite voltage difference V across the junction appears, 
and the phase difference starts to evolve in time as 

dt h 

1.2.1 Small Junction 

For a junction small enough to assume that the phase along it is constant, 
and in the limit of low temperatures (T — > OK), the two equations above yield 
quantum mechanical description of small Josephson Junctions (e.g. see |24j). 
In the quantum description, the phase difference A0 is replaced by a more 
general gauge invariant phase <fi conjugate with the number of Cooper pairs 
n transmitted across the junction, [(f), n] = i. The junction Hamiltonian is 
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Figure 1.2: a) rf-SQUID and its equivalent circuit, b) Its potential energy as 
a function of phase for two different ratios of inductive El and Josephson Ej 
energies (E L = &1/2L, $ e /$ = 0.5). 

determined by junctions' charging energy E c , and Josephson energy Ej: 



where C is the junction capacitance, n q is the externally induced continuous 
polarization charge, and $o is the magnetic flux quantum. 

Depending on the ratio between charging and Josephson energies the dy- 
namics of Josephson junctions can easily be analyzed in two respectfully ex- 
treme regimes. The limit Eq 3> Ej is the one of the Cooper pair box and it 
was discussed earlier. In the opposite limit, Ej ^> E c , the phase difference 
across the junction becomes a good quantum number that describes the quan- 
tized magnetic flux that is passing between the two electrodes of the junction. 
After the dropping of the constant Ej term from (jl.6|) . the Hamiltonian can 
be expressed in in the "coordinate" representation (n —> —id^) as 



Hjj = E c {n - n q f - Ej(l - cos(0)), 



(1.6) 




H = E c (id^ + n q ) 2 - Ej cos(0). 



(1.7) 
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The periodic potential enables expressing the solutions of this Hamiltonian in 
the form of Bloch functions 

where the continuous nature of n q gives rise to the existence of energy bands, 
while u(4>) has period of 2tt and can be expressed in the terms of Mathew 
Functions [25] . 

The use of the Cooper pair dynamics for design of charge qubits was obvious 
after one junction electrode was made into a superconducting dot that trapped 
Cooper pairs. The Josephson flux dynamics can also be used for design of 
qubits by trapping the flux that passes through the junction on one of its sides 
as proposed in [211] • This requires inserting the junction into a superconducting 
ring as shown in figure 11.2b . The flux that passes through the junction from 
the right gets trapped inside the superconducting ring and in this process it 
induces a persistent current in the ring. This way shorted junction is called 
rf- SQUID 2 . 

Putting a junction into a superconducting loop of inductance L gives rise 
to the inductive part in the Hamiltonian (jl.7j) . The loop shorts the junction 
and the induced charge n q becomes irrelevant since it can always be screened 
out. Generally, some external flux $ ext can be applied to the loop and thus in- 
ductively coupled with the Josephson phase. Then the rf-SQUID Hamiltonian 
is 

H = -E c d\ - Ej cos(0) + |^(0/2vr - <W$ ) 2 . 

The periodicity of ()1.7|) is clearly broken. Depending on the ratio of induc- 
tive energy El = $o/2L and the Josephson tunneling Ej, the shape of the 
potential varies from one with a single minimum to one with many local po- 
tential minima (figure fOb ). If El/Ej « 3 and $ ext /$o ~ 1/2 the potential 

2 SQUID is short for Superconducting QUantum Interference Device. The device is a 
good flux to voltage transducer and it is often used for precise flux measurements. These 
measurements require monitoring the device with tank circuit operating at radio frequencies, 
thus the prefix rf. 
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is characterized by two local minima (figure 11.3b ) each corresponding to a 
different flux state that can mutually tunnel through the barrier, and whose 
relative energy spacing can be controlled by outside applied flux. This way 
optimized rf-SQUID can be used as qubit. The experiments 3 j2Hl EH] have ob- 
served macroscopic quantum tunneling and the superposition of macroscopic 
persistent current states ffigure fOb h but the need for the large inductance L 
has hampered the observation of coherent oscillations since large inductance 
makes the qubit very sensitive to the fluctuations in the environment. As 
argued in [31], the need for the large inductance loop can be overcome by 
building a three-junction SQUID f figure H~3*h ) characterized by smaller exter- 
nal inductance and thus less sensitive to the outside fluctuations. The ability 
of the three-junction SQUID to preform as a flux qubit was demonstrated in 

1.2.2 Long Junctions 

If the junction interface is long, the passing of a fluxon between the electrodes 
cannot be considered instantaneous but rather as a real time event. In the 
certain limit of the junction parameters, the long junction supports the prop- 
agation of fluxons along the junction where the dynamics of the fluxons is 
affected by the injected or outside induced junction current. This way opti- 
mized long junction is often referred to as Josephson transmission line (JTL). 
At low temperatures these fluxons exhibits quantum dynamics j2J; an d they 
can be considered as individual quantum particles propagating in one or two- 
dimensional space, whose metric and the potential are fully determined by the 
junction parameters and the inserted current density respectively. 

The varying phase of a long junction can be analyzed by treating the 
long junction as a distributed structure. In the case of a long junction with a 
varying superconducting phase along only one if its interfaces, the junction can 

3 These experiments used rf-SQUID with the single junction replaced by dc-SQUID in 
order to control the tunneling in addition to the energy splitting of the rf-SQUID. For more 
on dc-SQUID see [231 E2 EE!- 
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Figure 1.3: a) Energy levels for rf-SQUID for E L /E C = 3, = 11/8tt, 

A/ El = 0.06. b) Experimental setup used in observing first superposition of 
quantum flux in a qubit and the plot of measured energy relative to the mean 
level energy as a function of bias e. The plot clearly shows energy splitting 
caused by the existence of tunneling. For more see [3U] c) Three-junction flux 
qubit and imbedded dc-SQUID magnetometer used for the first observation of 
flux qubit oscillations. 
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be subdivided into strips with the strips treated as small junctions connected in 
parallel by inductors as shown in the lower part of figure The Hamiltonian 
of this one-dimensional structure is given by 



where q n and n are conjugate phase and charge variables of strip n that satisfy 
[4>ni Qm] = 2ei8n,m, while C n , Ej n and Em are respectively the capacitance, 
Josephson tunneling, and inductive energy of the strips. 

The equation (jl.8|) can be expressed in the continuous approximation if the 
limit where the strip size is much smaller than the length of the junction but 
larger than the superconducting penetration depth and the spacing between 
the electrodes. The continuous Lagrangian density of the uniform junction 
with length expressed in the units of Josephson penetration depth Aj and time 
in the units of plasma frequency u> p , reduces to the sine-Gordon Lagrangian 



ej = %JJ2it, e L = (%/2n) 2 /l , u 2 p = ej/a, X 2 = e L /ej, 



where J c , cq and Iq are Josephson current density, junction capacitance density 
and inductance density respectively. 

In the framework of semi-classical quantization of sine-Gordon Lagrangian 
|33| . the propagating fluxon is represented as a relativistic soliton of size nXj. 
The average position and velocity of the soliton are specified by its two col- 
lective coordinates, while the other coordinates correspond to the internal 
quantum excitations of the soliton. At low temperatures (ksT <C ^eje L ), and 
as long as the velocity of fluxon is not on the order of Xjuj p , the fluxon cannot 
excite these internal degrees of freedom. It behaves as a relativistic quan- 
tum particle of mass M = yfleje^ / \ 2 jUJ 2 propagating free in one- dimensional 
relativistic space with "speed of light" given as c = Xjuj p . 






(1.9) 
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Figure 1.4: Long Josephson junction, depicted in the bottom of the drawing, 
can be viewed as distributed structure consisting of small Josephson junc- 
tions connected in parallel by inductors. In the certain limit of the junction 
parameters, the long junction supports propagation of fluxons in between its 
interfaces. It can perform as a detector, since the near-by flux qubit couples to 
it inductively and modulates the scattering matrix of the propagating fluxons. 

Coupling a flux qubit with the JTL (figure II .4j) creates potential in the 
JTL, and the Hamiltonian (jl.8|) and Lagrangian (jl.9j) get a source term whose 
explicit form depends on the type and strength of coupling. In any case, the 
source term is localized as long as the qubit size is much smaller than the JTL, 
and the freely propagating fluxons will now scatter off the potential caused by 
the localized source term and the transmission properties of the 1-D channel 
formed by the optimized JTL will depend on the qubit state. This setup can 
be used for determining the state of the qubit, but if not designed properly 
the propagating fluxons along the JTL will interact with the qubit and cause 
decoherence. 

1.3 Modelling Decoherence 

The qubits introduced above are affected by the changes in the environment. 
This is not just caused by the measurement of the qubit, rather it can be 
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contributed to many other factors. For example, 1/f charge fluctuations induce 
variations of the static Coopers in the pair box, while the inductive loop of 
the rf-SQUID easily couples the flux qubit with the outside electromagnetic 
fluctuations that vary the bias of the qubit. The coupling of the qubit with the 
outside noise sources adds additional degrees of freedom to the system whose 
random nature smear the coherent phase of the condensates and decohere the 
qubit. In many cases the microscopic origin of the decoherence is unknown. 
Therefore, to study the decoherence it is necessary to devise a model that treats 
the qubit as a quantum system in contact with some general environment which 
in the appropriate limit reduces to a stochastic process. 

System-plus-bath or system-plus-reservoir is the most natural way to ap- 
proach this problem. It considerers the weak entanglement between the system 
of interest S, and the much larger surroundings that act as a heat bath or a 
reservoir pMllHSIIHo] . Reducing the general equations of motion by tracing-out 
the environmental degrees of freedom leaves the equations dependent explic- 
itly only on S. In this process the environment is assumed classical, static, 
and unaffected by S. Thus any energy that moves from S is not stored, but 
lost in the environment. This results in slow, irreversible transfer of the en- 
ergy from S to the environment that is reflected in the reduced equations of 
motion. In the further works [371 EE1 EH] , the general reservoir was replaced 
by infinite number of harmonic oscillators. This took into an account that the 
environment by itself is a quantum system and relaxed the earlier condition 
of strictly classical environment. 

Spin-Boson^ is the common name for the model that precipitates after 
the system-plus-reservoir model is applied to a quantum two-level system |4Uj . 
Explicitly, it consists of the qubit Hamiltonian (jl.2|) linearly coupled to a bath 
consisted of a large number of harmonic oscillators 

h N N 

H SB = H+ -a z J2 Me* + 4) + Yl (1.10) 

i=l i=l 

4 The name originates from the analogy between two-level system and spin coupled to 
many oscillators that act as bosons. 
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where q and c\ are creation and annihilation operators, and Aj is the coupling 
between the i-th oscillator and the qubit. The environmental effects are phe- 
nomenologically represented by the noise spectral density S(u) of the coupling 



For weak coupling, the qubit will slowly couple to the oscillators in the 
bath, and in the large oscillator bath limit (N — > oo), the energy dissipated by 
the qubit will never come back. As a consequence, the qubit will relax to some 
energy equilibrium state determined by the qubit-bath coupling. This equilib- 
rium state will effectively become the ground state of the qubit and for this 
reason the qubit will become static and stop exhibiting coherent oscillations. 

1.4 Quantum Measurement 

In order to measure the qubit it is necessary to entangle it with a detector 
whose classical output signal is correlated with the state of the qubit. Because 
of the entanglement, a definite state of the output signal localizes the qubit in 
one of the states determined by the detector-qubit coupling. This is the same 
process as the reservoir induced decoherence discussed earlier. The major 
difference is that the extracted information from the qubit is not all lost, but 
that a part or preferably all of it is reflected in the detector output. 

Single-shot measurement is characterized by perfect correlation between 
the qubit and the detector output. This requires strong coupling between the 
detector and the qubit. Thus measuring the qubit extracts all the informa- 
tion and instantaneously collapses the qubit into the state correlated with the 
observed output. 

As an example we can consider the setup as shown in figure 11.41 where 
the JTL is situated close to the flux qubit consisting of appropriately opti- 
mized rf-SQUID. In the case of negligible decoherence and strong inductive 
coupling between the qubit and JTL, a fluxon propagating from the left along 




(1.11) 



i=l 



1(3 



the transmission line is totally reflected if the qubit is in \R) state, and totally 
transmitted if the qubit is in the \L) state. So after the scattering, the trans- 
mitted fluxon is entangled with the \L) and vice- versa. Detecting weather 
the fluxon was reflected of transmitted will collapse the qubit into the state 
entangled with the scattering outcome. 

In order to recover the statistical properties of the superposition state of 
the qubit, it is necessary to repeat the single-shot measurement on a certain 
number of the identically prepared and time evolved qubits The ob- 

servation of the quantum coherent oscillations of the qubit requires further 
measurement of the superpositions at different times. 

Coupling the qubit weakly with the detector does not result in instanta- 
neous wave function collapse, and for this reason the measurement provides 
only limited information about the qubit. This is reflected in the detector 
output which is not perfectly correlated with the qubit, rather it is charac- 
terized by existence of the overlap between the two signals that represent the 
corresponding qubit states. 

If the qubit and JTL in the earlier example are weakly coupled, the fluxon 
scattering matrix is only affected partially by the state of the qubit. The 
transmitted or reflected fluxon is entangled with both states of the qubit. 
Detecting the scattering outcome in this case does not collapse the qubit, but 
provides information in which state the qubit is more likely to be. 

The advantage of the weakly coupled detector is its ability to monitor 
qubit continuously for some time. During the measurement process, the qubit 
will slowly localize under the influence of the detector back-action which does 
not necessarily keep the qubit coherent, rather the off-diagonal density matrix 
element in the flux basis is limited from above as < y/ <J Lhit)® unit) ■ 

The back-action dephasing and the continuous extraction of information 
are two parallel dynamical processes where the rate of back-action induced 
dephasing is never smaller that the rate of information extraction. In the 
ideal case of the two rates being equal, the detector back-action keeps the 
qubit coherent, and the detector will extract information form the qubit as 
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fast as it is localizing it. Any deviation from the ideal measurement will 
result in loss of information since the decoherence rate will be larger than the 
information extraction rate. This is a general result that can be obtained from 
linear response theory j32] applied to a Hamiltonian H = H q + Hj + xf , where 
H q and Hd are general quantum system and detector respectively that are 
linearly connected through a general quantum system coordinate x and some 
detector "force" /. The detector provides output signal o(t) that consists of 
the detector nose q(t) and response Xx(t). 

Assuming that the detector is static, the dynamics of the measurement is 
determined by three correlators of detector parameters that describe detec- 
tor back-action (/(0)/(£)), detector response (f(0)q(t)), and the output nose 
(q(0)q(t)). In this framework, the condition of the ideal measurement is gen- 
eralized of having detector satisfy Re((q(0) f (t))) = 0. Furthermore, it can be 
shown that in the case of the ideal measurement of coherent oscillations the 
detector signal-to- noise ratio is limited to four from above [35] . 

This limitation is a direct consequence of the detector back-action and it 
can be traced back to the nature of the detector-qubit coupling. If x in the 
coupling is to be made a constant of motion of the measured system, then the 
detector would preform quantum non-demolition (QND) measurement that 
has no upper limit on the signal to noise ratio jH]. To illustrate QND mea- 
surement consider ffuxons propagating along the JTL detector weakly coupled 
to the flux qubit. If the fluxon insertion frequency equals to the qubit fre- 
quency, the passing fluxons will "kick" the qubit in exactly the same way and 
in this process exert no back-action on the qubit. This detector will thus be 
able to determine the flux state of the qubit, but it will also make the qubit 
static in the basis of detector coupling [35] . 
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Part I 
Quantum Decoherence 
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Chapter 2 



Weakly Coupled 
High-Frequency Noise 

The purpose of this chapter is to develop theoretical description of decoherence 
in the case of two coupled qubits subject to high frequency (high-f) noise. 
It follows a standard, Markovian approximation approach for description of 
weak decoherence based on the perturbation expansion of the density matrix 
evolution equation. The result obtained within this simple scheme is useful 
as a starting point for noise treatment of coupled qubit systems and as a 
benchmark for more elaborate models. 

2.1 The Model 

The model is a system-plus-bath model where an arbitrary number of quantum 
operators Qi of the quantum system H are each weakly connected with a 
different bath Hjg though a generalized bath force /j. The baths can be 
correlated. Furthermore, they are assumed to be much larger that the quantum 
system and initially in an equilibrium. The large size and the weak coupling 
makes the baths unaffected by the qubit. They remain in equilibrium and 
can be separated out from the total density matrix of the system. Explicitly, 
Ptot{t) = p(t)pB, where the two density matrices on the right respectively 
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describe the tensor product of quantum system and the total "bath-described" 
environment. 

If the system is initialized at t — 0, its Hamiltonian 



(2.1) 



expanded to the second order in the coupling yields differential equation that 
describes the evolution of the qubit's density matrix in the interaction picture 
as: 



In the equation (|2.2[) . the force operators can be factored out and grouped 
into their expectation values (fi(t)) and correlations (fi(t)fj(t')), where (■•■) = 
trsipB ' ' ' }• Since H can always be renormalized to absorb the non-diagonal 
force operators {(fi{t)} 7^ 0), it can be assumed that (fi(t)) = and that the 
environment is represented only by the force correlators. According to the 
central limit theorem (eg. see j3B]) and the the Wick theorem HO] using 
the second order correlators to describe the noise created by a large system is 
sufficient. The central limit theorem implies that the macroscopic size of the 
environment makes the noise generated by it Gaussian, while the higher order 
correlators of Gaussian noise are always expressible by Wick's theorem in the 
terms of second order correlators. Therefore, the correlators (fi{t)fj{t')) or 
equivalently their spectral densities 





f dt ' tr B{mmt), ifjWQiWipWM]} 





(2.3) 



fully describe the environment in our case. 
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The further assumption of Markovian 1 nature of the system, implies that 
p{t) depends only on p(t) and justifies replacing p(t') — > p{t) in the integral 
part of the equation (12.2(1 . The high frequency nature of the noise makes the 
bath response time instantaneous and the integration limit in (|2.2|) can be 
extended to infinity. The differential equation ()2.2|) simplifies to 

=-«E M{[Qi®,Q& -t')p(t)} </*(f )£(«)> 

-[Qi(*),P(*)Oi(*-O]</i(0)/i(O». (2-4) 

Expressing the operators Qi(t) in the interaction representation, (i.e. Qi(t) = 
exp(itH / h)Qi(0) exp(—itHo/h)) and then recognizing that all but non-oscillatory 
terms of the equation ()2.4j) average out, yields differential matrix equation that 
describes the relaxation of the quantum system. The equation is in general 
solvable but progressively more complicated with the increasing energy spec- 
trum of Hq. 



2.2 Single Qubit 

Long and extensive history of studying two-level quantum systems in optics 
and NMR calls for the discussion about the decoherence in a single qubit to 
be very short. The more detailed derivation of the results can be found in 
standard literature, i.e. [3H] . 

We start by specifying the two-level system coupled to a bath through its 
vertical polarization 

H = ~- [ea z + Aa x ) - -j-a z + H B , 

with the bias e, tunnelling A and noise / all having units of frequency. The 
choice of the coupling corresponds to the electrostatic interaction through finite 
coupling capacitance for charge qubits, or magnetic interaction for flux qubits. 
1 Evolution of a Markovian system does not depend on its history. 
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The interaction part of the Hamiltonian in the interaction basis is: 

Ht7 = M=^ \ea z - AaJ . 



2Ve 2 + A s 



After following all the steps outlined above, it is easy to arrive to standard 
two-level relaxation result that implies longest coherence times and shortest re- 
laxation times for qubits operating at zero-bias, (e = 0), point. The relaxation 
rate at this optimal point is T = 1/Ti = S(A)/2 while the decoherence rate is 
7 = 1/T 2 = 5*(A)/4. Consistent to Fermi's golden rule they are different by 
factor-of-2. 



2.3 Double Qubit 

Motivation for studying decoherence in coupled qubits is provided by the first 
successful double charge qubit experiment [20], where it was found that the 
decoherence rate for quantum coherent oscillations in two qubits at the optimal 
bias point is with good accuracy factor-of-four larger than the decoherence 
rate in effectively decoupled qubits. An interesting question for theory is 
whether this factor-of-four increase of the decoherence rate is a numerical 
coincidence, or it reflects some basic property of the decoherence mechanisms 
in charge qubits. As will become clear from the discussion below, the theory 
developed in this section favors "numerical coincidence" point of view. Other 
aspects of decoherence in coupled qubits have been studied before numerically 

in jH El Ell- 
in general, it is well understood that decoherence rates of different states 
of two coupled qubits can be quite different if the random forces created by 
the qubit environments responsible for decoherence are completely or partially 
correlated. Most importantly, in the case of complete correlation, the qubit 
system should have a "decoherence-free subspace" (DFS) spanned by the states 
|01), 1 10) [JEIEIIIS], since completely correlated external environments can 
not distinguish these states. In contrast, the subspace spanned by 1 00) and 
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1 11) experiences decoherence that is made stronger by the correlations between 
environmental forces acting on the two qubits. So the role of the quantitative 
theory in description of decoherence in the dynamics of coupled qubits is to 
see to what extent subspaces with different decoherence rates participate in 
the qubit oscillations for different methods of their excitation. 

2.3.1 The model and environmental correlations 

The Hamiltonian of the system of two qubits coupled directly by interaction 
between the basis-forming degrees of freedom is: 

H = hJ2 M° + *s°f + + MM 2) , (2-5) 

J'=l,2 

where ex's denote the Pauli matrices, v is the qubit interaction energy, Aj is 
the tunnel amplitude and Ej is the bias of the j-th qubit. Four energy levels 
of the Hamiltonian ()2.5|) are shown schematically in figure 12.11 as functions of 
the common bias e = E\ = £2 of the two qubits. This work considers quantum 
coherent oscillations in the qubits biased at the "co-resonance" point |2Uj . 
where £\= £2 = 0. Such bias conditions are optimal for the oscillations. 

It can be shown explicitly that the occupation probabilities of the qubit 
basis states (that are of interest for us) are insensitive to the phases of the 
qubit tunnel amplitudes Aj, so without the loss of generality we will assume 
that Aj's are real. The Hamiltonian (J2.5|) at the co- resonance reduces then to 

H = hJ2 ^ + • ( 2 -6) 

In the basis composed of eigenstates of the ax operators, the Hamiltonian 
(12. 6J) can be diagonalized easily. Eigenenergies and eigenstates are: 

e x = m, = |[(7 + v)(\oo) + in)) + 7)(|01) + |io»], 

E 2 = -Ml, |^ 2 ) = ±[(7 - 77)(|00> + |11)) + (7 + 77)(|01> + |10»], (2.7) 
E 3 = he, ^3) = |[(« + £)(|00> - |H>) + (a - /3){\10) - |01))], 
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Figure 2.1: Schematic structure of the energy levels of the two coupled qubits 
as functions of the common bias of the qubits. The inset shows the diagram 
of the decoherence-induced transitions between the levels at "co-resonance" 
point where the bias vanishes. 



E 4 = -he, |V> 4 > = - «)(|00) - |11)) + (a + /3)(|10) - |01»], 
where 

Q = (A 2 + v 2 ) 1 ' 2 , e = (6 2 + v 2 ) 1 ' 2 , 

and A = Ax + A 2 , 5 = A x — A 2 . The states \kl) with {A;, /} = {0, 1} in equa- 
tions (J2.7|) are the eigenstates of the operators a z in the natural notations: 
ajp\kl) = (-l) k \kl) and af ] \kl) = (-l) l \kl). 

We assume that external environments responsible for the decoherence cou- 
ple to the basis-forming degrees of freedom of the qubits. The interaction 
Hamiltonian is then: 

H z = hJ2m^ ] - (2-8) 

.7=1,2 

The random forces fi^it) acting on the qubits are in general correlated. To 
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describe the weakly dissipative dynamics of the system in the basis of states 
(J2.7I) induced by the interaction ()2.8|) with the reservoirs, we use the standard 
equation ()2.4j) . Proceeding in the usual way, we keep in equation ()2.4|) only the 
terms that do not oscillate, and therefore lead to changes in p that accumulate 
over time. Equations for the matrix elements p nm , n, m — 1, ...,4, of p in the 
basis ()2.7|) are transformed then as follows: 

Pnm = [ — Pnmi&mm ~ a nn)(^jj'(^) a rnm ~ ^jj'(^) a nn) 

-Pnm Yl^nl a kn^n'^ ~ € k) + a mk a ki^jA € m ~ e fe)) 
k 

+Snm Pkki^k^jf^k ~ 6n) + ^An ~ ^) 

k 

+ E PuJSo&iSjfa - e m ) + S}^ - e m ))) . (2.9) 
(k,i) 

Here a£L denote the matrix elements (n\ai^\m), the last sum is taken over 
the pairs (k, I) of states that satisfy the "resonance" condition: 

e k - ei = e n - e m , (k, I) ^ (n, m) , 

and 

poo 

Jo 

The first term in equation (|2.9|) represents "pure dephasing" that ex- 
ists when the system operators that couple it to the environment have non- 
vanishing average values in the eigenbasis. As one can see explicitly from 
equations (|2.7|) . the average values <Tnn of are vanishing in all states, so 
that there is no pure dephasing term in the evolution of the density matrix 
in our case. The fact that cxin are vanishing can be related to the vanish- 
ing slope of the system energies with respect to variations of the bias in the 
vicinity of the co-resonance point - see figure 12.11 Since all coefficients in 
the eigenf unctions (|2.7|) are real, the matrix elements <Jnm are also real. For 

(i) ~ 
real a nm -, imaginary parts of the noise correlators Sfj in the second term on 
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the right-hand-side of equation ()2.9|) represent the decoherence-induced shifts 
of the systems' energy levels. These shifts do not affect decoherence and we 
will neglect them in our discussion. With these simplifications, equation ()2.9j) 
takes the form 

Pnm = (2.10) 

[ - (Pnm/2) ^(^4°^ (e n - e k ) + ^a^ReS^m - e k )) 
i,j'=i,2 k 

+5nrn ^ Pkk^nk a kl ReS if ( e k ~ <*) + ^ Pkl^ V^S 'jf^ - € m )] , 
k (k,l) 

where 

/oo 
■oo 

The function S 12 characterizes the correlations between the environmen- 
tal forces acting on the two qubits. For instance, if the two qubits interact 
with different environments and vi, t>2 are uncorrelated, S12 = 0, whereas 
5*12 = S11 = S22, if the qubits are acted upon by the force produced by a sin- 
gle environment coupled equally to the two qubits. While the correlators Sn 
and S22 are strictly real, S12 can be imaginary, and = S12. Non- vanishing 
imaginary part of S12 corresponds to the non- vanishing commutator [/i,/2] 
and implies that the two qubits are coupled to the two non-commuting vari- 
ables of the same reservoir. While this is probably not very likely for qubits 
with the basis-forming degrees of freedom of the same nature (which in a typ- 
ical situation should be coupled to the same set of environmental degrees of 
freedom), the non- vanishing ImS'12 should be typical if the qubits have dif- 
ferent basis-forming variables, (e.g. see system proposed in [52]). Using the 
spectral decomposition of the correlators Sjj> (u) and Schwartz inequality, one 
can prove 2 that for arbitrary stationary reservoirs the correlators satisfy the 

2 Similar proof in slightly different context of linear quantum measurements can be found 
in P2]. 
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inequality that imposes the constraint on iSi 2 (a;): 



(2.12) 



If the reservoirs are in equilibrium at temperature T, the correlators also satisfy 
the standard detailed balance relations: 



Equation (J2.1U|) with the noise correlators (J2.11J) govern weakly dissipative 
time evolution of the two coupled qubits in a generic situation. Below we use 
them to determine decoherence properties of quantum coherent oscillations of 
the qubits. 

2.3.2 Quantum coherent oscillations in coupled qubits 

One of the most direct ways of excitation of quantum coherent oscillations in 
individual or coupled qubits that will be discussed in this work is based on 
the abrupt variation of the bias conditions [T%1 |2*U] . If the qubits are initially 
localized in one of their basis states, e.g. 1 00) , and abrupt variation of the bias 
brings them to the co- resonance point, the probabilities for the qubits to be 
in other basis states start oscillating with time. 

In the simplest detection scheme (realized, for instance, in experiment [2H] ) 
the probability for each qubit to be in the state |1) is measured independently 
of the state of the other qubit. Quantitatively, these probabilities pj are ob- 
tained from the projection operators Pf 



p^TripPj}, P 1 =Y,\lk)(kl\, P 2 =J2\ kl )( lk \- ( 2 - 14 ) 



Finding explicitly the matrix elements of Pj from the wave functions ()2.7|) . 




(2.13) 



fc=l,2 



fc=l,2 
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one gets: 



Pi(t) = - + {ar l + [3 1 )Re[e- i "- t (p 42 (t)-p 13 {t))) 

+ (a T -^)Re[e-^+*(p 14 (t)+p 32 (t))], (2.15) 

P2{t) = \ - (a 1 + (3 V )Re[e-^- t (p 13 (t)+p 42 (t))} 

+ (ar / -/? 7 )Re[e-^(p 14 (t)-p32(t))], (2.16) 

where u± = Q ± e, and as in the equation ()2.10|) . the matrix elements of 
the density matrix are taken in the interaction representation. Equations 
(|2.15|) . (|2.16p . and (|2.9|) show that the waveform of the coherent oscillations 
in coupled qubits is determined by the time evolution of the two pairs of the 
matrix elements of p: 

Pl3 = -ri 3 pi3 + M_p 42 , p 4 2 = -r 4 2P42 + «+Pl3 , ( 2 -17) 

p i4 = -r 14 p 14 + v_p 32 , p 32 = -r 32 p 3 2 + V + p U ■ (2.18) 

The decoherence rates in these equations are determined by the rates of 
transitions between different energy eigenstates: 

r 1 3 = ^(rl +) + r^ + r( +) + ri + )), 
r 14 = ^(rS^ + rS + ) + rr ) + r( +) ), 

r 42 = i(r« + r^ + rw + r«). (2.19) 

where labelling of the transitions is indicated in the inset in figure 12.11 Tran- 
sition rates are: 

rS ±) = Re^^,,(±^ + )4 ) 4' ) , 

3,3' 
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31 ; 



hi' 



42 5 



(i)^(j') 

32 



(2.20) 



3 J 



The superscripts ± refer here to transitions in the direction of decreasing (+) 
or increasing (-) energy. After finding matrix elements a nm from the wave 
functions (|2.7|) . we see explicitly that transitions between the states 1 and 2, 
as well as 3 and 4 are suppressed. Since the corresponding matrix elements 
are zero, the rates (|2.20j) are: 



r\ = 


— S\ 1 
2 


1 - 


8A + v 2 ' 
eft 


+ ^ 22 


1 + 


r 2 = 




1 + 


5A + u 2 ' 
eft 


+ \s 22 


1 - 


r 3 = 
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1 + 


5A + is 2 ~ 
eft 


+ \s 22 


1 - 


r 4 = 


-Su 


1 - 


5A + u 2 ' 
eft 


+ \s 22 


1 + 
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eft 
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eft 
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eft 
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eft - 



- Re5'i2 
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Re512I ftf ' (2 - 21) 
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The transfer "rates" u, v in equations ()2.17|) and (J2.18J) are: 

u±= E rEVS-M^) , w ±= E ^^^/(i^-) ■ (2-22) 



Explicitly: 



u = -S n 


1 - 


5A + v 2 ' 
eft 


" ^22 




1 4 


(5A + i/ 2 " 
eft 


+ \s 22 



5A-u 2 ' 

eft 
SA-v 2 

eft 



zlm5'i2- 
■ ilm5'i2 



VUJ_ 

ftT 



fte ' 



(2.23) 



Equations (|2.21jl and ()2.23)1 do not show the frequency dependence of noise 
correlators, which is the same, respectively, as in the equations ()2.20|) and 
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Each pair, ()2. 17j) and (J2.18|) . of coupled equations can be solved directly 
by diagonalization of the matrix of the evolution coefficients with a non- 
orthogonal transformation. In this way we obtain for the pair of equations 
(GUI): 

. , U 3 (0)(n + n_e^+* + c 2 e~ 7 -*) + cu_M0)(e~ 7+1 - e^-*)] 

Pl3(*) = ~ — ~ 2 , 

U + U- + c z 

M [p 42 (0)(u + u_e-^ + c 2 e-^) + c M+ p 13 (0)(e^+* - e"^)] 
P42{t) ~ + • (2 ' 24) 

where 

7± = (r 13 + r 42 )/2 ± [(r 13 - r 42 ) 2 /4 + 1/2 , 
c = (r 13 - r 42 )/2 - [(r 13 - r 42 ) 2 /4 + U+U J\ 1/2 , (2.25) 

and p 43 (0), p 42 (0) are the initial values of the density matrix elements that 
depend on preparation of the initial state. If, as in the experiment [20]; the 
qubits are abruptly driven to co-resonance maintaining the state |00), these 
initial values are: 

Pis(O) = ^(y + rj)(a + /3), p 32 (0) = J( 7 - 77) (a + (3) , 

p 14 (0) = i( 7 + r/ )(/3-«), p 42 ( ) = i( 7 -r / )(a-/?). (2.26) 

Another type of initial conditions that will be discussed in this work is starting 
the oscillations from the state |10). In this case: 

Mo) = \(j-v)(<x-P), fte(p) = \('y + v)(ot-P), 

p 14 (0) = ^(y-nXa + p), p 42 (0) = i( 7 + r/)(a + /5). (2.27) 

Equations ()2.22|) and (J2.23)) follow directly from the wavefunctions (|2.7|) : p nm (0) 
{n\i){i\m), where \i) is the initial state. 
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Solution of the other pair ()2.18|) of coupled equations is given by the same 
equations (J2.20|) and ()2.25|) with obvious substitutions: u± — > v±, r i3 — > T u , 
T42 —>■ ^32- In this work, we are mostly interested in the low-temperature 
regime ksT/h <C e, Q, when transitions up in energy can be neglected. In 
this regime, V- — > 0, and equations for the evolution of the density matrix 
elements are simplified. For instance, for u_ — > 0, c ~ u + u„/(Ti3 — T42), and 
equations ()2.24|) are reduced to: 

Pi 3 (t) = pi 3 (0)e- ri3t , 

p 42 (t) = p 42 (0)e- r ^ + r " + pi 3 (0)(e- ri3f -e- r ^), (2.28) 

1 13 — 1 42 

where now Ti 3 = (ri + T 2 + T 4 )/2 and r 42 = T 3 /2. 

Time evolution ()2.28|) of the density matrix elements together with the 
rates (J2.21|) and (J2.23j) . and initial conditions (|2.26|) and (|2.27|) determines the 
shape of the coherent oscillations in two coupled qubits. In the next section, 
we discuss this shape in several specific situations. 

2.3.3 Results and Conclusions 

The shape of coherent oscillations determined in the previous section depends 
on the large number of parameters: temperature, the degree of asymmetry of 
qubit tunnel energies and couplings to the environments, frequency dependence 
of the decoherence, and strength and nature of the correlations between the 
two reservoirs. We analyze some of these dependencies below. 

Experimentally-motivated case 

We start by considering the situation that is close to the experimentally real- 
ized case of oscillations in coupled charge qubits [20] • As argued above, the 
correlations between environments in this case should be real: ImSi2 = 0. 
The oscillations are excited by driving the system to co-resonance in the ini- 
tial state 1 00). Equations of the previous section give in this case the following 
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expression for the shape of the oscillations: 
Pi{t)= l -- l -[Ae^ t + B{e 

-IP 



-r 32 t 



+ D(e- Tl4t + 



Tl3 — 1^42 



Tl4 — 1^32 



r M t _ e ~ r42 *))] coscj_t 
e ~r 14 t _ g-rsat))] COSC( , + t ; (2.29) 



where 



5 A + z/ 2 z/ z/ 5 A + v 2 v v 

efi fie' efi fie' 

^ _ 5 A + v 2 v v D _ 5 A + v 2 v v 

efi fi e ' efi fi e 

Equation for P2{t) is the same with signs in front of S, u + and v + reversed. 
As a simplifying assumption we take Su = S 2 2 = S. The decoherence rates in 
equation (|2.29J) then are: 
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+ 2 5cK) 7k 



(2.30) 



SJu). 

2 v 



z/u; + 

l_ fiT' 



q/ ^A 5A 
u+ = -S[u+)— ,v+ = -S{u-)—, 

where S c (u) = KeSuiuj). To enable comparison of these rates to those of 
individual qubits, we note that the rate of decoherence of oscillations in a 
single qubit with vanishing bias is equal to S(Aj)/2 for the jth qubit. 

The functions Pi^if) for the qubit parameters, 5, v, and S, close to those 
in experiment j20] are plotted in figure l2~2l under additional assumption that 
the decoherence is the same at two frequencies, u + and cu_. 3 The curves are 
plotted for the two situations: when decoherence is completely uncorrelated 



3 The decoherence strength S was taken from the data for the single-qubit regime in |2()j . 
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Figure 2.2: Probabilities pj to find jth qubit in the state |1) in the process 
of quantum coherent oscillations starting with the state 1 00) of two coupled 
qubits. The decoherence strength is S — 0.08A. Solid and dashed lines cor- 
respond, respectively, to the decoherence that is uncorrelated (S c = 0) and 
completely correlated (S c = S) between the two qubits. 

(5*12 = 0) and completely correlated (Su = S) between the two qubits. One 

can see that the difference between the two regimes is not very big numerically. 

The correlations between the two reservoirs leade to the effective decoherence 

rate that is increased in comparison with the uncorrelated regime by roughly 

30-^50%, although the description with a total correlation is not quantitatively 

quite appropriate - see equations ()2.29|) and (|2.30|) . 

The increase of the effective decoherence rate by correlations illustrated in 

figure I2~21 can be related to the fact that the initial qubit state, |00), belongs to 

the subspace where the correlations increase the decoherence rate, despite the 

mixing of this subspace with the DFS 4 where the decoherence rate is decreased 

in the eigenstates (|2.7|) of the coupled qubit system. This implies that increase 

4 Here and below we use the term "DFS" for the subspace spanned by the 1 01) and 1 10) 
states, although for interacting qubits it, strictly speaking, does not fully have the properties 
of real DFS. 
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of decoherence rate by the correlations should be to a large extent insensitive 
to the qubit parameters. This conclusion is supported by the case of identical 
qubits (5 = 0), when u + = v + = and equation (|2.29|) is reduced to a very 
simple form: 

= \ ~ \[ i 1 + ^) e- Tlst cosuj.t + (l - £j e- r ^coscu + t] , (2.31) 

and P2{t) = Pi{t). One can see from equations (|2.30)) that both decoherence 
rates relevant for equation (|2.31|) . Ti3 and 1^32 increase with increasing cor- 
relation strength S c . Equation (j2.31|) shows also that the description of the 
oscillation decay with a single decoherence rate can be quite inaccurate: for 
weak interaction, v < Q the amplitudes of the two (high- and low-frequency) 
components of the oscillations are nearly the same while their decoherence 
rates can be very different. 

Excitation into the DFS 

Now lets discuss decoherence properties of the oscillations in coupled qubits 
in the case when they start with the initial qubit state |10). We note that in 
the case of experiment similar to j20], such an initial condition would require 
separate gate control of the two qubits, since the bias change bringing them 
into co-resonance is different in this state for the two qubits. Since the state 
1 10) belongs to the DFS in the case of completely correlated noise, one can 
expect that oscillations with these initial conditions will be more sensitive to 
the degree of inter-qubit decoherence correlations than oscillations with 1 00) 
initial condition, and that the effective decoherence rate will decrease with 
correlation strength. All this indeed can be seen from equations (|2.28|) with 
the initial conditions ()2.27|) that correspond to the 1 10) state. Under the same 
assumptions as were used in equation ()2.29|) . we get for the now different pi(t) 
and p 2 (t) : 

Pj (t) = \ - t±L{ [ Aje -™ + Bj (e-™ + -M^(e- r ^ _ e -r« 4))] CQSUJ _ t 

Z 1 1 3 — 1 42 
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+ [C 3 e- r ^ + D,(e- r ^ - i 1)3V+ (e" Fl4t - e~ r ^))] cos^ + t} , j = 1,2 (2.32) 



14 J- 32 
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Ai = 1 + 



d = l 

and the amplitudes A 2 , B 2 are given by the same expressions with 

5 -> -5. 

For identical qubits equation ()2.32|) reduces to: 

Pi(*) = ^ " (l + ^) e- r42t cos^t + (l - e^'cos^t] . (2.33) 

Which in the conjunction with equations ()2.30|) show that in contrast to equa- 
tion (|2.31|) . the decoherence rate of the low- frequency component that has 
larger amplitude is strongly suppressed by the non-vanishing inter-qubit noise 
correlations S c : r 42 = ~(l + v/£T)[S(u-) —S c (u-)]. This means that the shape 
of the coherent oscillations in coupled qubit starting with the state 1 10) should 
indeed be more sensitive to the strength of these correlations than the shape 
of the oscillations starting with the 1 00) state. The conclusion also remains 
valid in the case of not fully symmetric qubits as one can see from figure 12.31 
which shows the shape (J2.32|) of the 1 10) oscillations for the same set of ex- 
perimentally realized parameters as in figure l2~2l Even in this case, there is a 
pronounced weakly decaying component of the oscillations if the decoherence 
is completely correlated between the two qubits. For partial correlations, the 
effective decoherence rate is reduced. 

In summary, we have developed quantitative description of weakly dis- 
sipative dynamics of two coupled qubits based on the standard Markovian 
evolution equation for the density matrix. This description shows that deco- 
herence properties of currently realized oscillations in coupled qubits are not 
very sensitive to inter-qubit correlations of decoherence, while relatively simple 
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Figure 2.3: Probabilities pj to find j th qubit in the state |1) in the process 
of quantum coherent oscillations starting with the state 1 10} of two coupled 
qubits. Qubit parameters are the same as in figure 12.21 Solid, dotted, and 
dashed lines correspond, respectively, to the decoherence that is uncorrelated 
(S c = 0), partially (S c = 0.5S), and completely (S c = S) correlated between 
the two qubits. 
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modification of the excitation scheme for the oscillations should make them 
sensitive to these correlations. 

Before moving on, lets briefly discuss the applicability of our approach to 
the realistic Josephson-junction qubits. As we saw above, one of the main 
features of equation (|2.10|) is that the pure dephasing terms disappear at the 
co-resonance point and the remaining decoherence is related to the transitions 
between the energy eigenstates. Therefore, within the approach based on 
equation (|2.1U|) . the decoherence rates are on the order of half of the transition 
rates. On the contrary, the experiments with charge qubits (see, e.g., [T§] ) 
indicate that decoherence rates are larger than the transition rates even at 
the optimum bias point when the pure-dephasing terms should disappear. 
Apparently, this is related to the low-frequency charge noise |53J UH] that is 
coupled to qubit strongly enough for the lowest-order perturbation theory in 
coupling ()2.4j) to be insufficient. This implies that the theory presented so far 
might be only qualitatively correct for realistic charge qubits, and that a more 
accurate non-perturbative description of the low-frequency noise is needed in 
order to achieve quantitative agreement with experiments. 
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Chapter 3 



Non-Perturbative 
Low-Frequency Noise 

The experimentally observed decay time T 2 of coherent oscillations is typi- 
cally shorter than the energy relaxation time Ti even at optimal bias points 
[T§1 EH EH where the perturbation theory predicts factor-of-two differ- 
ence between the two times and no pure dephasing terms. Furthermore, the 
observed increase in two-qubit decoherence rates (201 cannot be explained by 
the perturbation theory results of the previous chapter. Qualitatively, the ba- 
sic reason for discrepancy between 7\ and T2 is the low-frequency noise that 
can reduce T2 without changing significantly the relaxation rates. Mecha- 
nisms of low-frequency, or specifically 1/ / noise exist in all solid-state qubits: 
background charge fluctuations for charge-based qubits |H3 ESI IS3 EE] , and 
impurity spins or trapped fluxes for magnetic qubits jSHl- Manifestations of 
this noise are observed in the echo-type experiments jHSl- Low-frequency noise 
for qubits is also created by the electromagnetic fluctuations in filtered control 
lines. 

The aim of this chapter is to develop quantitative theory of low-frequency 
decoherence by studying single and double qubit dynamics under the influence 
of Gaussian, noise with small characteristic amplitude Vq, and long correlation 
time r. In the case of single qubit we will show that the expression describing 
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the decay times of the coherent qubit oscillations is a non-perturbative result 
whose strength is controlled by the zero-frequency noise spectral density 2vqT. 
For long correlation times r ^> A -1 , where A is the qubit tunnelling amplitude, 
2-UqT can be large even for weak noise Our analytical results are exact 

as function of 2v%t in this limit. 

In the second part of this chapter, the same non-perturbative technique 
is applied to double qubit system operating at the co- resonance point. It is 
shown there that the decoherence rates of the coupled qubits compared to 
the decoherence rates of the individual uncoupled qubits qualitatively double. 
However, the change in the energy spectrum in going from uncoupled to cou- 
pled system plays an important role that quantitatively varies the qualitative 
factor-of-two change of the decoherence rates in much wider range. 

3.1 Single Qubit Decoherence by Low-f Noise 

We start with the standard qubit Hamiltonian with the fluctuating bias v(t), 

H = ~[Aa x + (e + v(t))a z ], (3.1) 

where the noise v(t) has characteristic correlation time r. Therefore, its cor- 
relation function and spectral density can be taken as 

2v; 2 r 

(v(t)v(t')) = vle-^ , S v {u) = , (3.2) 

1 -+- (LJT) 

where Vq is the typical noise amplitude in units of radial frequency and ( ■ ■ • ) 

denotes average over different realizations of the noise. We assume that the 

temperature T of the noise-producing environment is large on the scale of the 

cut-off frequency 1/r, and can be treated as classical. 1 

The decoherence is a decay of the off-diagonal element of the qubit 's density 

matrix in the interaction basis. Therefore, evaluating decoherence is equivalent 

1 In the regime of interest, 1/t«A, the temperature can obviously be still small on the 
qubit energy scale. 
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to evaluating time evolution of the a operator in the same basis. In the path 
integral representation this can be expressed as: 



where ^[o -- ] represents integration over all possible paths a~ can take in going 
form <J~(0) to cr~(t), angled brackets represent noise averaging, and T and T* 
are respectively forward and reverse time ordering operators. 

One contribution to the expression above is from the transitions between 
the two energy eigenstates. The transitions are caused by the high-frequency 
part of the noise spectrum and they can be described by the means of the 
perturbation theory as shown in the earlier chapter. The condition of weak 
noise Vq -C A makes the transition rate small compared both to A and 1/r 
ensuring that the perturbation theory is sufficient for the description of transi- 
tions. The additional effect of weak noise dynamics affecting the qubit (|3.1|) is 
"pure" or adiabatic dephasing. As discussed qualitatively earlier, the fact that 
the noise correlation time is long, r ^> A -1 , makes the perturbation theory 
inadequate for the description of pure dephasing. For low-frequency noise, a 
proper (non-perturbative in Vqt) description is obtained by looking at the ac- 
cumulation of the noise-induced phase between the two instantaneous energy 
eigenstates in the expression ()3.3|) . 

Since the transitions correspond to averaging over the all possible states of 
cr~, they can be omitted from the path integral ()3.3|) by neglecting the integra- 
tion over cr _ . This simplification reduces the functional integral over a com- 
plicated Keldysh contour to one that is considerably simpler (fig. I3.1JI . After 
putting the explicit values for the Hamiltonian H = — ^/A 2 + (e + v(t)) 2 a z /2, 
and noting that a z and a~ anti-commute, we can define the factor F(t) that 
describes the time-dependent, low frequency suppression of coherence between 
the two states as 




(3.3) 




(3.4) 
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Figure 3.1: The general path integral of (a~(t)) has a complicated average 
that connects forward and backward Keldyish time contours (left). Neglecting 
the transitions, simplifies the integral to two independent and equal averages 
over each of the time directions (right). 



If vq <C A, one can determine the rate of accumulation of this phase by 
expanding the energies up to the second order in noise amplitude v(t), i.e. 



'A 2 + (e + v(t)) 2 = tt+ (ev(t'))/n + (A 2 v 2 (t'))/(2tt 3 ), tl 2 = A 2 + e 2 , and 
evaluating ()3.4j) in the interaction representation. The pure dephasing term 
is: 

m =\ exp vL [-¥ + ^ 1 \ dt })- (3 - 5) 

For Gaussian noise, the correlation function (|3.2|) determines the noise 
statistics completely. In this case, it is convenient to take the average in 
equation (|3.5|) by writing it as a functional integral over noise. For this purpose, 
we start with the "transition" probability p(t>i, t> 2 , <5t) [BUI EI] for the noise to 
have the value v 2 time St after it had the value V\. 

, Art 1 / (V2 ~ Vye-V') 2 1 nR , 

p( "" 1,2 '' i)= ^,g(i- e -^) exp r 2^ ( i- e -^) ) • < 3 - 6 > 

Using this expression we introduce the probability of specific noise realization 
aspo(v 1 )-p(v 1 ,V2,5t 1 )-p(v2,v 3 ,5t 2 )-..., where p (v) = (2ttvq)~ 1/2 exp{-w 2 /2uo} 
is the stationary Gaussian probability distribution of v, and SU are some small 
discrete time steps. Putting the explicit values for p's given above with the 
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same 5t into the formula for F(t) we obtain the following expression 

exp{-vf/2v 2 } 



F(t) 



(27rvg)( JV+1 )/ 2 (l - e -2«/r)-Ar/ 2 

N 

JJexp 

n=l 



X 



[v n+1 - v n e 



-6t/r\2 



2v 2 (l - e~ 28t /-] 



exp < — i5t 



2„,2 "1 



A 2 w 



n 2Q 3 



After expanding up to second non-vanishing order of 5t and taking the limit 
5t — > 0, the above expression can be written as F(t) = S(X)/S(1), where S(X) 
is the following functional integral: 



S(X) = (const) x - / dv(0)dv(t)Vv(t f ) 



K 



v(o) 2 + v(t) 2 r* , .. 2 2 n . , 

x exp <; — — / dt(v 2 + v 2 + 2iav) 

K Jo 



with k, a and z given as 



(3.7) 



K= 1 + 



2u; 2 A 2 r\ 5 
fi 3 



a 



r 



(3.« 



Since the average in equation ()3.5|) with the weight (|3.7J) is now given by the 
Gaussian integral, it can be calculated in the usual way (i.e. see |62j): 



F(t) 



F (t)exp 



-a 2 ( — - 2[coth — + d" 1 
r 2t 



e*/ 2r [cosh(Kt/r) + l±^! s inh(^/r)]- 1/2 . 



(3.9) 



3.1.1 Single Qubit Results and Conclusion 

Equation ()3.9|) is the main analytical result for dephasing by the Gaussian 
noise. To analyze its implications, we start with the zero-bias, (i.e. e = 
0) case, where pure qubit dephasing vanishes in the standard perturbation 
theory. Dephasing ()3.9|) is still non-vanishing and its strength depends on 
the noise spectral density at zero frequency S v (0) = 2v\r expressed through 
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k = a/1 + is, s = S v (0)/A for convenience. The numerical simulations show 
that equation (J3.9|) and the transition contributions fully account for qubit 
decoherence as shown in figure 13.21 For small and large times t equation ()3.9jl 
simplifies to: 



F(t) 



1+t/i 



1 +t/r + ist/2r 
2V«e- (7+ ^7(l + 



lV2 



(3.10) 



where 



7 



1 

27 



5 2Nl/2 



1/2 



(3.11; 



Besides suppressing the coherence, the noise also shifts the frequency of qubit 
oscillations. The corresponding frequency renormalization is well defined for 
t > r: 



,2)1/2 _ ! 



lV2 



Suppression of coherence (|3.10|) for t <C r can be qualitatively understood 
as the result of averaging over the static distribution of noise v. In contrast, at 
large times t ^> r, the noise appears to be (^-correlated, the fact that naturally 
leads to the exponential decay (|3.10j) . This interpretation means that the two 
regimes of decay should be generic to different models of the low-frequency 
noise. Indeed, they exist for the non-Gaussian noise considered jHSlEll; and 
are also found for Gaussian noise with 1/ / spectrum . Crossover between 
the two regimes takes place at t ~ r, and the absolute value of F(t) in the 
crossover region can be estimated as (1 + s 2 ) -1 / 4 , i.e. s determines the amount 
of coherence left to decay exponentially. The rate (J3.11j) of exponential decay 
shows a transition from the quadratic to square-root behavior as a function 
of S v (0) that can be seen in figure 1331 The figure also shows that the decay 
rate extracted from numerical simulations of Gaussian noise agree well with 
the theoretical predictions for quite large noise amplitude v . 

Non-zero qubit bias e leads to the additional dephasing F(t)/F (t) de- 
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Figure 3.2: The envelope of total decoherence due to adiabatic dephasing and 
the transitions as predicted by equation ()3.9|) and the golden rule (dotted line) 
plotted together with Monte Carlo simulation of the quantum oscillations for 
the same qubit (solid line). The simulation consists of 10 5 realizations over 
noise parameters specified on the plot. 
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Figure 3.3: The rate 7 of exponential qubit decoherence at long times t ^> r 
for e = and noise with characteristic amplitude Vq and correlation time r. 
Solid line gives analytical results: Eq. (J3.11|) . Symbols show 7 extracted 
from Monte Carlo simulations of qubit dynamics. Going from left to right 
along the horizontal axes the number of Monte Carlo realizations gradually 
changes from 10 7 to 10 4 . Inset shows schematic diagram of qubit basis states 
fluctuating under the influence of noise v{t). 
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scribed by the last exponential factor in equation (|3.9|) . The contribution 
from Fo(t) is of the same form as in zero-bias case but now with the replace- 
ment s — > s(A/Q) 3 . The additional dephasing exhibits the crossover at t ~ r 
from "inhomogeneous broadening" (averaging over the static distribution of 
the noise v) to exponential decay at t ^> r. In contrast to zero-bias result, the 
short-time decay is now Gaussian: 



In 







L^b(*)_ 


Q 2 1 



w 2 rt/(l +is(A/fi) 3 ), i»r. 



We see again that the rate of exponential decay depends non-trivially on the 
noise spectral density S v (0), changing from direct to inverse proportionality 
to S v (0) at small and large s, respectively. 

The approach outlined above can be used to calculate the rate of exponen- 
tial decay at large times t for Gaussian noise with arbitrary spectral density 
S v (u)). Such noise can be represented as a sum of the noises (|3.2jl and ap- 
propriate transformation of the variables in this sum enables one to write the 
average over the noise as a functional integral similar to (J3.7j) . For calculation 
of the relaxation rate at large t, the boundary terms in the integral (J3.7j) can 
be neglected. The integral is then dominated by the contribution from the 
"bulk" which can be conveniently written in terms of the Fourier components 



(2/t) 1/2 / dt'v(t') sin u n t' , u n = im/t, 



and, (•••) = / Dv ■ ■ -exp{ — (1/2) J2 n \ v n\ 2 / S v {u n )}. Combining this equation 
and equation (|3.5|) we get at large t: 



F{t) = exp 



e 2 QS v (0) 1 
_W + iS v {0)A 2 + n 



du In ( 1 + i 



fi 3 



.(3.12) 



For unbiased qubit, e — 0, this equation coincides with the one obtained by a 
more involved diagrammatic perturbation theory in quadratic coupling [Bo] . 



47 



3.2 Double Qubit 



Now we proceed to determine the adiabatic dephasing decoherence rate for a 
coupled qubit system operating at the co-resonance point. As described earlier 
the model double-qubit Hamiltonian at co-resonance point is: 

ffo = fiE &A i} + toWo® + v ti)^\ (3-13) 

3=1,2 3=1,2 

where again cr's denote the Pauli matrices, v is the qubit interaction energy, 
Aj are the tunnelling amplitudes of the two qubits, and Vj are two weakly 
coupled, in general correlated noise sources. As it can be seen from figure 
12.11 applying the idea of pure dephasing being dependent on noise induced 
variation of the energy splitting between the levels of interest, is simplified at 
the co-resonance point by the vanishing linear dependence between the energy 
levels and the noise parameters V\$- 

After following the same arguments as for the single qubit case of neglecting 
the transition rates and only concentrating on the variation of energy split- 
ting between the corresponding energy levels, the general expression for pure 
dephasing expressed a path integral in the interaction picture is: 



Pij(t) 




dt'v m c^lv n ) (3.14) 



PiM 

(ij) E {(13), (14), (32), (42)}, 

where c^fm are elements of 2 x 2 symmetric matrix describing the dependence of 
the energy splitting as a function of noise sources obtained by Taylor expansion 

ciij) _ 1 PjEj-Ej) 

n > m 2h dv m dv n ' 1 1 

The averaging in ()3.14|) is now done over two noise sources whose dependence 
on integration variable t' is implied. 

The analytic values of the coefficients (J3.15)) can be easily obtained with 
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the help of the perturbation theory. As expected from the vanishing first- 
order energy variations at the co-resonance point, the first-order corrections 
to energy vanish, i.e. 

&\ = o. 

The second order deviation in energy is non-vanishing, and it can be shown 
that 

d 2 E l _ 2 y ma l j i) \^)(^ j \ai m) m 

dv n dv m 4^ Ei~ E i 

where the wave functions are defined in (|2.7j) . The coefficients of four 
matrices are explicitly given in table 13.11 
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c 22 
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eQ(Q - e) 
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eQ(Q + e) 
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eQ(Q + e) 


eQ(Q + e) 
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Table 3.1: The coefficients of symmetric 2x2 matrix c^' for four different 
density matrix elements pij. 

Lastly, we are left to discus the complications that arise in evaluation of the 
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dEj 
dv lt z 



dH 

dvto 



expression (J3.14J1 due to the correlation between the noise sources V\ and V2- 
In general, correlation between any two stochastic processes can symmetrically 
be expressed through a coefficient A such that 

vx(t) = _(vT+Ai?i(t) + vT^i? 2 (t)), (3.16) 
v 2 (t) = -j= (yi + Xvxit) - vi-a^w) , 

where v\ y2 are two uncorrelated noise sources, i.e. {v\(t)v2^S)) = 0. Therefore, 
after an appropriate orthogonal transformation ()3.16|) . any averaging over two 
correlated noise sources can be expressed as averaging over two uncorrelated 
sources. 

Assuming that the two noise sources &12 are identical, then (vi(t)vi) = 
(^2(^)^2) — Vq€~^^ t , and the final expression for adiabatic dephasing of the 
density matrix elements needed to calculate the two probability values pi,2(t) 
defined in JQH) is: 



Pij(t) 



exp (1 * m Cg£v n \) (3.17) 



m,n=l,2 

(ij) E {(13), (14), (32), (42)}, 



Vl,2 
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= (ciY } 


+ ^-2^(1- 
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/nf(ij) 
°12 


_ /nr(ii) 
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= (c& ) -cSVl- 


-A 2 



Each noise source is fully described by conditional probability ([3.6)1 . and the 
averaging (|3.17|) is straight forward along the lines of single qubit case. The 
un-normalized action is: 



S(C (ij) ) = (const) x 



50 



/ JJ dv m (0)dv m (t)Vv m exp<- 

m=l,2 I m=l,I 



4«8 



(3.18) 



exp 



4 ^o jo 



E 

m=l,2 



f; 2 

c 2 +^ 

m 1 r 2 



2vk; 2 



m,n=l,2 



An orthogonal transformation, u 1)2 — * t> 12 that diagonalizes the matrix 
decouples (I3.18j) into two quadratic integrals identical to those of optimally 
biased qubit (|3.9jl . where are in this case the functions of - i-e. the 
eigenvalues of matrix C^'. The constant term in each of the two integrals 
over the paths T>[v] get divided out by the normalizing terms consisting of 
interaction free action (J3.18|) . i.e. S'(l2x2)- All put together, the adiabatic 
contribution to dephasing of double qubit system operating at co-resonant 
point is: 
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F (t,K™)F (t,4 
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The expressions from the table yield the following eigenvalues 
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3.2.1 Double Qubit Results and Conclusions 

As we see from the decoherence rates of the double qubit contain two 

contributions from the low frequency noise that are identical to the one ob- 
tained for the single qubit (|3.9|) . Therefore, the conclusions of the single qubit 
carry over for the double qubit case. Specifically, there are the two regions of 
the decay - one when t r that is described by averaging over the statical 
distribution of the noise, and the other, i >t, where the noise appears to be 
delta-correlated and that is characterized by exponential decay. The expres- 
sions of the two decay coefficients at opposite time limits are given by (J3.1Uj) 
and (J3.11j) . Numerical simulations of the coupled qubit dynamics show very 
good agreement with the theoretical predictions of the decoherence rates as 
shown in the figure l3~4l 

Although the double qubit system contains twice as many terms that sup- 
press its coherent oscillations, we have to be careful not to understate the 
importance the rearrangement of the energy spectrum plays in the relative 
values of these coefficients and naively assert that the double system decoher- 
ence rate doubles as compared to the system's individual qubit decoherence 
rates. This dependence is rather non-trivial even under all the assumptions 
about the system and the noise considered so far. 

To illustrate this point we will consider the experiment [20] • Under the 
assumption of uncorrelated, identical baths, the low-f dephasing terms will be 
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tv/2n 

Figure 3.4: Log plot of numerical Monte Carlo simulations (solid lines) and 
theoretical decay predictions IY#(dashed lines) for evolution of |pi 3 | density 
matrix element for a coupled qubit operating at co-resonance point and with 
two noises exhibiting different levels of correlation (A = 0,0.5, and 1). The 
decay rates Tth are sum of high-f contributions from equations (|2.20j) and 
(12.21)) and low-f contributions (|3.19|) in the limit t > r as given by (|3.11|) . 
For the plots A = 5z/, 5 = 3u, r = 100/V, and initial condition 1^(0)) = |00), 
while the MC simulations consist of 10 5 realizations. 
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the same for qubit density matrix elements (13) and (42) as well as for the 
terms (32) and (14). At this point, it is advantageous to define a coefficient 
which represent the ratio of the large-time, (t ^> r), decoherence rates of the 
coupled qubit density matrix elements pij to the decoherence rates of the two 
individual, uncoupled qubits k — 1,2: 



m _ 7(S(Q)cf") + 7 (S(0)4 



* w = " 5 '-wwt) ' (3 ' 20) 



n ' 2r 



h ^v 2 + i\ 1/2 _ 1 



2 

(ij) 6 {(13), (14), (23), (24)}, k = l,2, 

where are given in ()3.19|) and A^ are individual qubit tunneling ampli- 
tudes. The the eight ratios ()3.20|) with the parameters given in [20] are plotted 
in figure l3~5l as functions of zero- frequency noise spectral density, S(0), in the 
units of A = Ai + A 2 . The factor-of-four increase in decoherence that was 
observed in the experiment can be explained by the "numerical coincidence" 
of zero-frequency noise spectral density being in the range that enhances the 
decoherence rates of the coupled qubits. 

The exact value of ratios (|3.2(J|) for the experiment [20] can be confirmed by 
determining S(0). Since the tunneling amplitudes of the qubits are different, 
the value of 5(0) within this limit of the model can be extracted from any two, 
out of two Ti and two T 2 times. Furthermore, if only the times of the same type 
are known, i.e. both T\ or both T 2 , they need to be measured precise enough 
so that their difference shows outside the margin of error. Unfortunately [2*U] 
provides only the two T 2 times that are the same within the margin of error 
and the exact ratio (|3.2(J|) cannot be determined. 

In summary, the non-perturbative treatment of the now-f noise from a 
fluctuators coupled to the basis-forming degrees of freedom has provided the 
analytical expressions which can fully within this model account for the in- 
crease of the decoherence rate beyond the predictions of Fermi's Golden Rule. 
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Figure 3.5: The ratio of low-f decoherence rates with the parameters taken 
from experiment and plotted as function of S(0)/A, where A = Aj + A 2 . 
The non-trivial relation between the change in the ratio of the decoherence 
rates between coupled and individual qubits in a double qubit Hamiltonian is 
evident. The expected factor-of-two increase of the decoherence rates of the 
coupled qubit as compared to the single qubit is not realized. Rather, it can 
vary over a much larger range. 
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In general, this is attributed to the accumulation of the phase in the qubit 
due to adiabatic low-f variations of the bias. Furthermore, the factor-of-four 
increase of decoherence rate of coupled qubits as compared to the individual 
qubits [2U] can be attributed to changes of the energy spectrum of the system 
when the two qubits are coupled as supposed when they are not. This supports 
"numerical coincidence" explanation to the decoherence rate increase rather 
than the existence of new decoherence mechanism. All the results have been 
verified numerically. 
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Chapter 4 



Realistic Quantum Modelling of 
Systems Subject to General 
Noise 

The ideas of adiabatic dephasing depending on the variation of the energy level 
spacing as the function of coupled weak noise and the transition rates being 
influenced by the value of spectral density of the noise at resonant frequency 
can be extended to a general n-qubit case in order to determine the desired 
decoherence rates of the system. Never-the-less, obtaining analytical results 
forces often a general system to be constrained in many ways such as operating 
at particular Hamiltonian constant in time and subject to specific, analytic 
noise spectral density 

The realistic applications often call for insights in to the dynamics of the 
system with time-varying parameters and subject to multiple noise sources 
described by some general spectral densities. In order to arrive to the result 
under these quite general conditions, it is necessary to abandon analytical 
expressions that were useful in qualitative understanding of the noise and 
resort to numerical simulations. The statistical nature of noise makes Monte 
Carlo (MC) methods ideal in realizing these simulations. 
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4.1 Simulating the Basic Qubit Dynamics 



In the simulation of the qubit evolution all the time-dependent qubit (e, A) 
and noise (v, u) parameters of the system Hamiltonian 

H(t) = -\ [(e(t) + v(t))a z + (A(t) + u(t))a x ] , (4.1) 

are considered constant during a given time step St, and the evolution of the 
qubit 's density matrix p can be calculated exactly during the time interval: 

p{t n + St) = exp(-iStH n /h)p(t n ) exp(iStHl/h), (4.2) 

where H n = H(t n ). At the boundary, the final values of density matrix are 
used as initial values of density matrix for the next interval with up-to-date 
parameters, 

p(t n+1 , H n+1 ) = p(t n + St, H n ). (4.3) 

Repeating this process eventually evolves an initially specified density matrix 
of a qubit to some final value along the path specified by deterministic Hamil- 
tonian and a random realization of each of the noises. Restarting this process 
over many different realizations of the noises and averaging the density matrix 
of the system at each time point yields complete time-evolved MC-evaluated 
density matrix of the system. 

The uncertainty of the averaged p(t) can be attributed to three sources: 
(1) computer truncation, (2) the assumption of constant parameters at the 
discrete time steps and, (3) statistical uncertainty due to the Monte Carlo 
averaging. 

Since the density matrix parameters are all on the order of unity the trun- 
cation error is at least sixteen orders of magnitude smaller. Its accumulation 
does not pose any serious problems since it is not compounded over the dif- 
ferent realizations and at the extreme one realization has 10 10 calculations. 
Thus even in the worst case, the result will contain error that is six orders of 
magnitude smaller and therefore negligible. 
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The error due to the discretization of the qubit parameters during the 
interval St is more serious because of the possibility of having unstable time 
step St. This would lead to numerical result diverging away from the exact, 
unknown value. The most direct way to improve on this error is to reduce 
St or use a mid-point or some more elaborate extrapolating method when 
converting the continuous parameters to discrete ones. In the case of unstable 
time step St, the simulation of the system without noise would produce density 
matrix exhibiting non-unitary behavior. For this reason, a very good check 
for an adequate choice of St consists of evolving the system without the noise 
forward in time from p(0) to p(tfi na i), then reversing the time parameter, (i.e. 
St — > —St), and evolving p(tfi na i) backward in time to p(0). The unitarity 
of quantum mechanics implies that the two matrices p(0) and p(0) should be 
identical if St is chosen appropriately. 

With St chosen such that p(0) = p(0) the only error in the result comes 
form the uncertainty due to Monte Carlo averaging. This error is given as 

Af J (jpUtn))) ~ ((Pi^W ,,,, 
&{Pij(tn)) = y jf ' ( 4 - 4 ) 

where N is the number of realizations, and the double-angled brackets denote 
arithmetic average over the realizations JU^J . Taking into account that the ele- 
ments of density matrix are on the order of unity, the error due MC averaging 
can be made reasonably small by increasing the number of realizations. 

The non-trivial nature of generating the noise form a given spectral density 
of its correlator does not allow the improvement to the variance ()4.4|) by uti- 
lization of stratified or importance sampling techniques. The 1 / \fN reduction 
of the variance can be improved by using the antithetic technique or by em- 
ploying quasi-random sampling. The antithetic technique [HZj is implemented 
by using a single realization of the noise twice - once as generated by random 
precess and second time reused but with opposite sign. The quasi-random 
sampling assures better covering of the sample space by generating the ran- 
dom numbers with the help of number theory that assures the "filling" of the 
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empty space between already sampled "random" points. This method can 
reduce the 1/y/N dependence tol/N [So] . 

4.2 Generating Noise 

Unlike in the case of previous qubits, the numeric Hamiltonian (|4.1j) has noises 
affecting both the tunneling and the bias of the qubit. This is often the case 
in the realistic situations. For instance in the experiment j3U] the tunneling 
noise can be attributed to the impurities in the layer separating the super- 
conducting electrodes of the Josephson junctions or to the variations in the 
modulating flux <& x d. c . as defined in figure 11.31 In the same experiment the 
bias noise originates form the variations in the applied flux <& x and from the 
dc-SQUID magnetometer. With exception to the impurity caused noise, the 
noises caused by the sources above are in general specified by some noise spec- 
tral densities Si(u) obtained form the known impedances of the circuits Yi(u) 
and the fluctuation-dissipation theorem [TT]. Thus, the random noise in the 
qubit simulation must be generated in a way that it matches the given spectral 
density. Since the significant portion of noise comes from the detector, which 
is usually coupled to the basis forming degree(s) of freedom of the qubit, the 
discussion that follows assumes that the only noise affecting qubit is coupled 
through z component so that u = in (|4.5[) . 

As discussed in jBH], the random variable v (t m ) generated from some given 
S(lu) as 

M 

v(mAt) = ^ y/AS(n5u)5u cos(n5umAt + </> n ), m = 0..M, (4.5) 

n=0 

is a Gaussian, and in the limit N — ► oo it has a correlator whose spectral 
density is S(u). In the equation (|4.5jh 8uj = l/(MAi) is discrete frequency in- 
terval while 4> n are uniformly distributed random phases taken randomly from 
the interval [0, 2ir). In the numerical implementation of the same equation, the 
Fourier transform is replaced by fast Fourier transform (FFT). In doing so, it 
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is important to choose M and At so that the cutoff time (MAt) is larger than 
the desired simulation time, and the cutoff frequency uo c exceeds the qubits 
maximal transition frequency. It is easily shown that u c = 1/ At [HE] . Further- 
more, in order to avoid numerical aliasing caused by approximating Fourier 
transform by FFT it is necessary to set the cutoff to be factor-of-two 1 larger 
than the maximal level spacing of the qubit VL max j^Hl EH] . 

The discrete Fourier transform given in ()4.5|) can be expressed as a real 
part of more general complex transform 



where the random phase is now included in the harmonic coefficient. This way 
written transform is more favorable for FFT computation because the random 
phase can be factored out. Also, real and imaginary parts of noise z(mAt) are 
uncorrelated since cos(nx) and sin(mx) are two orthogonal functions. Pre- 
forming one complex FFT yields two realizations of the noise, thus speeding 
up the simulation significantly. The antithetic variance reduction technique 
discussed earlier is used at this point in which case there are total of four noise 
realizations after single fourier transform (j4.fij) . 

Lastly the noise realization given in the time steps of At needs to be ex- 
pressed in the steps of the simulation St. This can be achieved along the same 
lines of discretization of the continuous parameters of the Hamiltonian (14.51) 
discussed at the beginning of this chapter. 

The approach outlined so far is very powerful, and it can produce results 

impossible to produce by analytical methods. Hadamard transform of a qubit 

state is a such example, and it's simulation shown in the figure l4~Tl 

lr This is not a strict limit. It was arrived to, through various computer simulations as 
well with the help from the references jHBJ 02] ■ 



M 




m = 0..M, (4.6) 



n=0 
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Figure 4.1: Hadamard transformation modelled by varying the bias as depicted 
in the graph a). The plots b) and c) show the variation of the respective density 
elements without noise while the plots d) and e) show the same evolution 
subject to noise. The simulation assumes A to be constant, while S{io) was 
supplied numerically. 
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4.3 Inclusion of Temperature 



Up to this point the temperature has not entered into the discussion, so the 
noise affecting the qubit is high-temperature noise with ksT/h 3> A, e. This 
implies that noise-induced up and down transition rates are the same. Realisti- 
cally this is not the case since qubit experiments are performed on temperature 
scales of the order of lOmK <C A, e, where the up rate is severely suppressed. 

As shown in the second chapter, the transitions are easily treated by per- 
turbation theory and they depend on the high-frequency part of the spectral 
density. Thus truncating the spectral density at the frequencies larger than 
some frequency u tr that is on the order and also below minimal qubit level 
spacing Vt min (figure I4.2j) eliminates the transitions from the simulation as 
shown in figure 14.31 It is important to realize here that although the values of 
spectral density above the truncation frequency are zero, the value At = 1/uj c 
in (|4.6|) stays fixed. If this was not the case, the removal of the transition fre- 
quencies form the generated noise will not be complete because it will contain 
aliased contribution. Aliasing is a serious problem in the signal processing 
|69j . The best way to assure that the aliasing is negligible while at the same 
time keeping in mind that increasing cut-off is prolonging the computation 
time, is to pick up lowest cut-off frequency of the noise spectrum for which the 
truncated spectral density produces no energy relaxation, i.e. T\ — > oo. 

Since the analytic expressions for the decay rates that obey the detailed 
balance relation are known (e.g. see [70]), the temperature dependent transi- 
tions are accounted for after adding- in "by hand" the appropriate coefficients 
in the equation ()4.2|) . Working explicitly in the interaction basis of the qubit, 
the diagonal density matrix elements acquire "by hand" inserted decay that 
restores the transitions removed by the truncation of the spectral density. The 
same can be said about the high-f decay component of the off-diagonal matrix 
element. Unlike before, the restored transition rates can be made to obey the 
detailed balance relation with an arbitrary temperature. The explicit modifi- 
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Figure 4.2: In order to eliminate transitions from the simulation the original 
spectral density a) needs to be truncated at the higher frequencies uj tr < n min 
b). 



cations to the equation ()4.2|) are: 



P00(*n + St) 
P0l(tn + St) 

The rates are given as 
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where T(Q) is effective temperature defined through FDT as 



2tt 2k B T(u) 



(4.7) 
(4.8) 



(4.9) 



(4.10) 



The inclusion of temperature requires additional specification of effective tem- 
perature T(u) or environmental resistance Re(Y(u)) over the transition fre- 
quency range. The identical simulations run at various temperatures exhibit 
increasing decoherence rates with the increase of the effective temperature. 
The figure 14.41 shows temperature dependence of Rabi oscillations of excited 
qubit state caused by the rf-pumping (insert). The direct relation between 
temperature and the degradation of coherent oscillations is obvious. 
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Figure 4.3: Two simulations with bias a) e = and b) e — 2A each run with 
original (solid line) and truncated (dashed line) spectral densities. The lack 
of the relaxation, a phenomenon solely caused by transitions, is apparent for 
truncated spectral density, while the relaxation is present in the simulation 
with original spectral density. 




Figure 4.4: Plot of Rabi oscillations at the optimal point (e = 0) simulated 
with different temperature values and the same rf- varied bias (insert). 
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4.4 Generalization to Larger Systems of Qubits 



The method outlined above for simulating the dynamics of qubit can be ex- 
tended to more complex qubit systems under the influence of classical noise. 
Inclusion of the temperature dependent transitions like as done above is im- 
possible since the perturbative analytic expression for time decay of the sys- 
tem's density matrix elements do not exist. In addition, the dimension of the 
Hilbert space scales exponentially (2^) to the number of the qubits N, and 
the density matrix approach is increasingly inefficient due to need for multiple, 
large-matrix multiplications. These shortcomings could be overcome by use 
of stochastic wave function approach [7TJ E21 EH EH EH] that is unfortunately 
beyond this work. 
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Part II 
Quantum Measurement 
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Chapter 5 



Continuous Weak Measurement 
with Mesoscopic Detectors 

5.1 Quantum Measurement and Mesoscopic De- 
tectors 

Textbook examples of quantum measurements often recall the Stern-Gerlach 
experiment where the measured phenomenon is quantum mechanical in na- 
ture, but the measurement itself is purely classical and could only determine 
the amplitudes of the system's wave function if repeated on a large number of 
identically prepared systems. This is a direct consequence of the design of the 
measuring equipment where the measured system, i.e. the beam of the elec- 
trons, is absorbed by the detector - a photographic plate or florescent screen. 
For this reason, the quantum state is always destroyed by the measurement. 

These destructive measurements are not general representative of the quan- 
tum measurement process as formulated by von-Neuman and the others [7H1 
EZ1 EHl EH ED] • The postulate of quantum measurement asserts that the mea- 
surement of an observable quantity q which gives a value qj leaves the density 
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matrix of the measured system p in the state 



1 



\Qi){Qj\p\Qi){Q: 



(5.1) 



Tr {\<lj)(<lj\p} 



where \qj) is the wave function of the corresponding measured value qj and 
the tracing in the denominator is done over both systems. Thus, measuring a 
general quantum state of a system and obtaining some result qj reduces the 
density matrix of the measured system in a way that the obtained result is now 
incorporated into the total, updated density matrix, where \qj) are defined by 
the detector-system coupling. 

In practice, the detector needs to be capable of measuring the state of the 
system while being at the same time minimally detrimental to the measured 
system's coherence. Since the coherence is purely quantum phenomenon, its 
preservation imposes necessity of measuring the quantum system by entangling 
it with another system that also exhibits some form of dynamics on the same 
energy scale. Furthermore, the need of having the result in the form of a 
classical, measurable signal makes mesoscopic phenomenon such as Coulomb 
blockade or mesoscopic transport in single or multichannel conductors prime 
candidates for a such task. 

By entangling a such mesoscopic detector with the measured system, the 
observation of a detector outcome does not need to destroy the measured sys- 
tem. If the detector is perfect, the reduction of the measured system's matrix 
is given by (|5.1|) . In general, the detector can act back onto the quantum 
system and the general state after the measurement can be represented as 



where £l(qj) describes the total back action of the detector on the measured 
system. If £l(qj) = \Qj)(Qj\ density matrix of the measured system will remain 
pure and the detector will be preforming in ideal or quantum-limited regime. 
The exact form of the operators Q(qj) depends on the type of the detector. 
If the measurement outcomes are orthogonal, i.e. \q%){qj\ = <% then \qj) 
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(5.2) 



Tr{\<lj){<li\Psvs} 
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are eigenvectors of the measurement basis and the measurement preformed by 
the detector will be projective or von-Neumann measurement. Assuming the 
quantum limited regime for now, the state of the quantum system after the 
measurement of qj will be projected into state given in (J5.1|) . As a consequence 
of reduction of the density matrix and the orthogonality of \qi)(qj\, any further 
measurements of the same quantity will give the identical outcome. 

Orthogonality condition implies the direct correlation between the quan- 
tum system and the detector. Achieving this requires strong coupling be- 
tween the two, which in general induces strong detector back-action beyond 
the limit of the ideal measurement. Reducing the back-action of the detector 
on the quantum system can be achieved by making the detector-system cou- 
pling weak. This leads to imperfect correlation between the detector output 
and the measured observable which is reflected in non-orthogonality of the 
operators \qi)(qj\- Relaxing the orthogonality constraint also implies that two 
successive measurements of the static quantum system do not need to give 
the same outcome, since two or more different outcomes can be implied by 
the same signal, i.e. the states \qj) are now some superpositions of the mea- 
surement basis states. In this case the output does not tell what state was 
measured, but it reflects the likelihood of the system being in the states that 
all could with some weights belong to the observed output. The successive 
weak measurements of the quantum system collapse its wave function to one 
of the states determined by the system-detector coupling. The slow reduction 
of the wave function enables continuous measurement of the quantum system 
and avoids the possibility of abrupt energy relaxation that is characterized by 
projective measurements which can result in total destruction of the quantum 
state due to the secondary reasons [HT] . 

A way to study weak measurement is in the framework of linear response 

theory where the interaction of weak detector force / and the quantum system 

coordinate x is assumed to be linear. 1 The Hamiltonian is specified as H = 

1 In certain situations the lowest order detector-system coupling is not linear but the 
analysis is carried along the similar pcrturbative steps outlined below, e.g. see |82| . 
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H q + Hd + xf, which in addition to the interaction coupling contains individual 
Hamiltonians of the quantum system H q and the mesoscopic detector Hd- 
Large size of the detector justifies the assumptions that the detector is static 
and that the changes in the detector are observable in the form of some classical 
output signal o(t). To the first order correction, the detector output can 
be written as consisting of the detector's noise v(t) and measurement signal 
superimposed over it, 



Since different times correspond to different and independent measurements, 
the noise of the detector is white and is completely specified by its correlator 
{v(t)v(t'))d = 2ixS v 5{t — £'), where S v is the noise spectral density that is 
assumed constant at the relevant low-frequency range of the output signal. 
The delta function is a consequence of the instantaneous detector response. 
The detector acts on the measured system through the detector force f(t) 
which induces stochastic motion to the quantum system. In this light, the 
correlator of the force describes the back-action of the detector, and it as 
well is characterized by some spectral density Sf constant in the relevant 
frequency range, it responds instantaneously and is independent of its values 
at the different times, thus {f(t)f(t'))d = 2nSf8{t — £'). The average detector 
force component is assumed to be zero, {f)d = 0, since any non-vanishing 
value can be absorbed by renormalizing the detector Hamiltonian Hd- The 
commutator of the cross-correlator between the noise and the detector force 
z ([/(£'), v(t)])d = h\S(t — t' — 0) essentially describes the information extraction 
from the measured system by the detector. The infinitesimal shift in the 
detector response is required in order to impose the causality in ([5.3)1 . 

In the case of the continuous measurement, the extraction of information 
from the measured system is not instantaneous but it is spread out over some 
time r m . Complementary to the extraction of the information from the quan- 
tum system the detector's back-action is dephasing the measured system's 
wave-function. In the case of static quantum system (H q = 0), the density 




(5.3) 
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matrix p of the measured system evolves in the eigenbasis of the observable q 
as: 



Pif = —fi(qj-Qj')f (* ) pjj' > ( 5 - 4 ) 

whose solution is given by Pjf(t) = Pjj'(0) exp(—i(qj — qj/)/h f f(t')dt'). After 
the averaging over the different realizations of f(t) and using the well known 
result, (exp(x)) = exp(x 2 /2), the decay of the density matrix elements due to 
the detector back-action is: 

Pn'it) = p jjf (0) exp(-r d i), T d = (5 . 5) 

At the same time, the dc component of the detector output signal changes 
by the value 5o = \{qj—qf) where A is the detector response coefficient defined 
earlier. In order to read-out the output signal, the background noise needs to 
subside to at least the half value of the 5o. Since the average background 
detector noise Ao averaged over some time interval At decays inversely to it 
as Ao = (2tzS v / At) 1 / 2 , the measurement time r m is: 

8nS v S f 

rm "[Ate-^))r { } 

From (|5.5jl and (|5.fij) the measurement time and back-action dephasing rate 
are related as 

r m T d = 8(7i/h\) 2 S q S f . (5.7) 

As detailed in ^2], the detector force and noise spectral densities can be 
defined as vector products since they are assumed constant at low frequencies: 



S l (u) = ~[S u (u) + S u {-«>)] 
S lk = (l\k), l,ke{v,f}. 
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Then the response coefficient can be expressed as A = 47iIm(Sf v )/h and the 
Schwartz inequality S vv Sff > \Sf v \ 2 , yields that 

lAI^^^^-.Re^) 2 ] 172 . (5.8) 

Combining (|5.7jl and (|5.8jl establishes a fundamental relation between the mea- 
surement time and back-action decoherence rate, 

T m r d > -. (5.9) 

Equation (|5.9|) implies that the acquisition of information is never faster than 
the dephasing due to back-action. At best, the quantum system can be mea- 
sured as fast as it is being dephased. This occurs for the detectors for which 
Re[Sf v (u)] = 0, where explicitly Sf v = f drexpiiuit) {fv{r)) d- The vanishing 
real part of spectral density function between the noise and the back-action 
implies that the detector noise is not correlated with the detector force thus the 
detector back-action is solely caused by exchange of information and not by 
the detector excitations of the measured system. This-way optimized detector 
is known as quantum-limited or optimal detector. 



5.2 Quantum Point Contact as a Detector 

The discussion above can be illustrated with a specific example of quantum 
point contact (QPC) detector situated in the proximity of the Cooper pair 
box qubit. QPC is a one dimensional channel whose resistive properties in 
the tunneling regime are very sensitive to electric fields [SHI EH EE EHl IHTj . 
Therefore, the current through biased QPC is modulated by charge present in 
the box. 

If the variations in the transmission of QPC are linear, then the detector 
and the interaction Hamiltonians are respectively: 
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H d = J2H E ^ ( 5 - 10 ) 

j=l,r i 

H = U= •/'''/• 

j,j'=l,r 

where a\ r and a^ r are electron creation and annihilation operators of the left 
and the right QPC electrodes, and are energies of electrons populating ith 
level. The coupling U is caused by the changes in the electron scattering poten- 
tial ±U (x) of QPC, where ± depends on weather the cooper pair occupies the 
box or not. Since the qubit is in a stationary state, its Hamiltonian is equiv- 
alent to zero. Left and right wave function of the electrons in the electrodes, 
ty(x)i, r , define the scattering matrix elements as Uij = J ^/*(x)U(x)^j(x)dx. 

Under the assumption of the QPC bias V larger than the qubit parameters 
and the environment temperature (eV » he, HA, fc^T), the average current 
passing through QPC is classical and in the zeroth order the correlators are 

(U(t)I(t + r)) = ^(t5T + u)5(r-0) (5.12) 
(I{t)I(t + r)) = 5(t) + ^<^(*)>, (5.13) 

where T and R are transmission and reflection probabilities, 5T is the change 
in transmission due to location of the cooper pair, while u is a dimensionless 
parameter that does not affect the current through the contact. Rather, it 
reflects asymmetry of the coupling between the box and QPC |4*3*] . 

The back-action dephasing can be obtained by solving an equation of mo- 
tion for the reduced qubit density matrix. Since the qubit is static, only its 
off-diagonal elements are of interest. In this case, the transition rate through 
the point contact depends on the qubit state and after the averaging over U 
as done in going form (|5.5|) to (|5.6|h the off-diagonal density matrix element 
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of the qubit evolves according to 

eV(ST) 2 + u 2 

Poi = -r>oi, T d =- . (5.14) 

The asymmetry factor u does not affect the current response, therefore it has 
an effect of increasing decoherence rate without any increase of the data ac- 
quisition rate T m = (SI) 2 /4:So, where Sq = 2e 3 VTR/nh is shot-noise spectral 
density and 51 = e 2 5TV/hn^\ EE]- The ratio of two decoherence rates, 

T m (ST) 2 

TT 1 = TT^TT 2 < 1, 5.15 

again is such that information extraction is always limited from the above by 
the qubit decoherence. 

Since the asymmetry factor does not affect the current response, it has no 
effect on information extraction from the qubit. On the other hand, the qubit 
decoherence (|5.14|) is increased by non-zero u and the ratio ()5.15|) is the most 
optimal for u — 0. This illustrates the interplay of back-action decoherence 
and extraction of information - only ideally optimized QPC (u = 0) preforms 
quantum measurement that extracts information from the qubit as fast as 
it is dephasing the quantum state of the qubit by its back-action. For non- 
optimal case (« ^ 0), a part of the information about the qubit is contained in 
the phase differences between the forward and backward scattering electrons 
in the QPC. Since this information cannot be extracted by measuring the 
ballistic properties of the QPC but by rather some "elaborate", or better 
said un-achievable interference between the forward and backward scattering 
channels, the information is lost and the QPC detector will dephase the qubit 
before it is able to measure its state completely. It is important to point out 
that the limitation (J5.15J) is independent of the strength of the qubit-detector 
coupling, rather it is an universal limitation of quantum mechanics. 
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5.3 Conditional Measurement 



In the regime of the quantum-limited detection, the overall evolution of the 
detector and the measured system is quantum-coherent and the only source 
of the information loss is averaging over the detector. For a detector, different 
outcomes of the evolution are classically distinguishable, and it is meaning- 
ful to ask how the measured system evolves for a given detector output. In 
the quantum-limited regime, specifying definite detector output eliminates all 
losses of information, and as a result there is no back-action dephasing present 
in the dynamics of the measured system conditioned on specific detector out- 
put. 

Conditional description in the quantitative form is obtained (see, e.g., |881 
I^IH]) by separating in the total wave function the terms that correspond to 
a specific classical outcome of measurement and renormalizing this part of the 
wave function so that it corresponds to the total probability of 1. 

In the conditional measurement approach the measured system is specified 
by its wave function |\&) = Conditional that the single detector 

particle was transmitted/reflected the amplitudes at are updated as 




(5.16) 



where tj,Tj are the values of the detector transmission/reflection as functions 
of the state j^) of the measured system. 

It is important to stress that the changes in the coefficients a, for a system 
with vanishing Hamiltonian (as we assumed from the very beginning) is un- 
usual from the point of view of Schrodinger equation, and provides quantitative 
expression of reduction of the wave function in the measurement process. 

If the dynamics of the detector is not quantum limited, then the informa- 
tion is lost and the dephasing is non-vanishing even in conditional evolution. 
To generalize equations (|5.16|) to the case of finite dephasing, we need to look 
at the change of the density matrix of the measured system due to a single 
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detector event. If the detector particle is specified as some free wave superpo- 
sition = J dkb(k)\k), then after tracing out the detector degrees of freedom 
from the total, system-plus-detector density matrix, the density matrix pij of 
the measured system changes as 

Pa - [ dk\b(k)\%(k)t*(k))/J2PffT f (5.17) 

f 

Pij - jdk\b(k)\\r l (ky j (k))/Y,PffRf (5-18) 



where 



are total transmission and reflection probabilities for a given state \qj) of 
the measured system. From these relations it is straight forward to see that 
quantum-limited detection requires the density matrix elements to remain pure 
in the conditional evolution. 



5.4 Measurement of Quantum Coherent Os- 
cillations 

In order to study the measurement of the quantum coherent oscillations, the 
output signal correlator (o(t + r)o(t)) needs to be evaluated explicitly for a 
measured system in non-static case. In the specific example of QPC detector, 
the correlator is the current- current correlator given in (J5.13|) . The first part 
of the expression is the contribution form the shot noise while the second part 
is the detector response to the qubit oscillations. To evaluate this {a z a z {r)) 
correlator, it is necessary to solve the equation of motion of the reduced density 
matrix of the qubit, which after substitution p = 1/2 + a yields: 

doo = AIm[<7 i], doi = (ze - T d )a 01 - iA<j n , (5.19) 
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where the back-action decoherence rate is given in (J5.14J1 and the normalization 
condition an = — croo was used to simplify the second equation. This equation 
is in general solvable, and in the case of zero bias (e = 0), the output signal 
spectral density S(u) = 2 J (/(0)/(r))e^ T rfr is: 

TA 2 (5I) 2 

where Sq = 2e 3 VTR/7rH is the shot noise spectral density. 

The zero-bias point is the most optimal for the qubit oscillations. For this 
value, the spectral density is peaked at the qubit oscillation frequency A. The 
maximal value of the peak above the shot noise is S max = (5I) 2 /T, and S(uj) 
is limited from above such that signal-to-noise ratio is 

Smax _ 4(5T) 2 ^ ^ ^ 



S (5T) 2 + u 2 

This again illustrates the limitation between the information extraction and 
the back-action decoherence. In the static case the limit was exhibited through 
the relation between the rate of decoherence and the rate of information ex- 
traction, while in the dynamic case it is exhibited through the signal-to-nose 
ratio that is limited to four from above. As before, the limitation is fundamen- 
tal (i.e. independent on the qubit-detector coupling). It represents the fact 
that in the time domain the oscillations are drowned in the detector's shot 
noise. 

The analytic expression for the output spectral density is possible after 
solving (J5.19j) with Cardain's formula, but the cumbersome expressions are 
not very insightful. Numerical evaluations of the spectral density for finite 
bias qubit can be found in 
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5.5 Quantum Non-Demolition Measurement 



Even if ideally optimized, the quantum measurements discussed so far do not 
exceed the standard quantum limit given by (J5.15j) or (J5.21j) . This is a conse- 
quence of the general choice for the system-detector coupling. For this reason, 
the measurement signal contains information about both conjugate variables, 
[p, q] = ±ih of the quantum system. The Heisenberg uncertainty about them 
(ApAq > h/4), is contained in the output signal and it prevents the signal 
from reflecting one of them exactly. 

The need to measure the mesoscopic observables beyond the quantum limit 
has produced a special group of measuring schemes known as quantum non- 
demolition (QND) measurements [""""1 19*01 l""Ti I92| . They all are based on the 
idea of measuring only one, i.e. q, of the conjugate variables of the system. 
This way constructed detector still acts back and perturbs the unobserved 
conjugate variable p, but since the detector is only coupled to q the output 
signal will not exhibit any consequences of the back-action on p. In other words 
the detector is coupled to the constant of motion of the measured system. 
Mathematically this implies that the measured observable and the operator 
Cl(qi) defined in (|5.2j) commute [75*j . 




Specifically, in the qubit measurement discussed so far the signal-to-noise 
limit (J5.21|) of monitoring quantum coherent qubit oscillations demonstrates 
the inability of the measurement scheme to obtain the information about the 
qubit beyond the quantum limit. The origin of this limit can be traced to the 
nature of the qubit-detector coupling that is only sensitive to a z projection of 
the qubit in some general, non-diagonal basis. For this reason, the detector is 
measuring directly the oscillating coordinate and consequently it localizes it. 
This is reflected through appearance of additional "quantum-limiting dephas- 
ing" which constrains the signal-to-noise spectral density of the output signal 
to four from above. Two known ways of overcoming this obstacle and realizing 




[?,ft(fc)] = 0. 



(5.22) 
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the QND measurement of qubit dynamics revolve around an idea of making 
the qubit-detector coupling time-dependent |H| I45j . 

The first of the designs solves this problem on the most straight-forward 
way - by coupling the detector to the unbiased qubit such that the detector 
follows the qubit as it oscillates in z^-plane. If the coupling is oscillating with 
some frequency uj c the system is described by Hamiltonian: 

1 / 

H = H d - ~ 2 [cos(u c t)a z + sm(u c t)cr y } . (5.23) 

In the regime of small detuning 7 between the measurement frame oscillation 
uj c and the oscillation of the qubit A, linear response language yields the 
spectral density of the output signal to be [H] : 

at , X2 T * + 1 - fa ~ A) 2 

S{u) = S + - j— 7—— , 7 = p , 5-24 

2huj 2 + (r e + 7J 2 r 

where A is the linear response coefficient, T e is the environment induced sup- 
pression of the coherence, and V is the detector back-action decoherence rate. 
In the case of zero detuning, 7 = and zero environment decoherence, T e = 0, 
the signal to noise ratio diverges for dc values of the output signal and the 
limit of standard quantum measurement (|5.15J) is surpassed. 

A somewhat different scheme introduced in jlH] achieves the QND mea- 
surement by employing the kicked-qubit technique where the charge qubit is 
coupled to QPC with the periodically switched coupling. In this case the de- 
tector has access to the qubit only at approximately discrete times. Ideally, the 
coupling part of the total Hamiltonian is a sequence of periodic delta functions 
separated by r and of strength //2, 

1 00 f(t) 

H = ~- [ea z + Aa x ] + £ 5(t - nr)^-a z + H d . (5.25) 

n=— 00 

The density matrix of initially pure quantum state \^) = a\R) + (3\L) after n 
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kicks with r = 2ir/VL is 



Pip) 



\a\ 2 Za{3* 
Z*a*p \(3\ 2 



) 



(5.26) 



where Z = exp[in(f) — n((5f) 2 )/2] is decay factor obtained by central limit 
theorem [2S|- It is obvious that the measurement scheme is preserving the 
information about the a z projection of the qubit, while the information about 
the remaining two projections a x , a y is lost. The spectral density of the detec- 
tor output is discrete, but again it diverges at dc frequency value. 

If the "kicked" QND measurement technique is applied to unbiased qubit, 
(H q = Aa x /2), then it is possible to synchronize the measurement pulses so 
that the detector interacts with the qubit only when the off-diagonal elements 
are identical to zero. Ideally in this case, the measurement will not affect the 
qubit evolution at all. Realistically, the coupling is not instantaneous, and the 
detector interacts with the qubit when the off-diagonal elements are nearly 
zero. In this process, the detector induces some small relaxation still beyond 
quantum limit decay rate that eventually collapses the qubit. More about this 
will be said in the following chapter. 
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Chapter 6 



Rapid Single Flux Detector 



Many flux qubits employ the measurement scheme based on the ability of a 
flux qubit to modulate the rate of tunneling out of the stationary supercurrent 
carrying states of a Josephson junction or a SQUID. This process can be viewed 
as tunneling of a magnetic flux quantum and has several attractive features 
as the basis for measurement. Most important one is the sufficiently large 
sensitivity which comes from the strong dependence of the tunneling amplitude 
on the parameters of the tunneling potential controlled by the qubit. There is 
an important disadvantage to this approach. In simple few-junction systems, 
the supercurrent decay brings the system into finite- voltage state characterized 
by large energy dissipation. This strongly perturbs the system and the detector 
itself, and makes it impossible to repeat the measurement sufficiently quickly. 
This strong perturbation destroys the state of the measured quantum system 
by changing it on an uncontrolled way. Both of these factors prevent the 
realization of the non-trivial quantum measurement strategies introduced in 
the previous section. 

The goal of this chapter is to suggest and analyze the flux detector based on 
the tunnelling of individual magnetic flux quanta, which avoids the transition 
into the dissipative state after fluxon-qubit entanglement by utilizing their bal- 
listic motion 1 . The detector has quantum-limited back-action on the qubit and 
1 Another way of avoiding the transition into the dissipative state in the course of mea- 
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time resolution sufficient for realization of non-trivial quantum measurements 
of superconducting qubits. 

6.1 Josephson Transmission Line as a Flux De- 
tector 

The detector is based on the ballistic motion of the fluxons in the Josephson 
transmission line (JTL) formed by unshunted junctions with critical currents 
I c , and capacitances C coupled by inductances L as shown in figure IdTI The 
detector can be viewed as the flux analog of QPC detector introduced in the 
previous chapter. Both detectors are based on the ability of qubit to control 
the ballistic motion of the independent particles (electrons in QPC and fluxons 
in JTL detector) through one- dimensional channel - the role played by JTL in 
the case of JTL detector as discussed in the introduction. Since the injection 
of fluxons can be controlled, the detector can be optimized to perform new 
measurement schemes (e.g. time-delay measurement) otherwise impossible 
with QPC. 

In the flux detector, the flux (f)^ e \x) generated by the measured system 
creates potential U(x) for the fluxons moving along JTL (figure lUTIj) . The 
potential U(x) is localized in some area on the JTL which acts as a scatter- 
ing region for the incident fluxons. The fluxons are one-by-one, periodically 
with period 1//, injected by the generator and their scattering characteristics 
(transmission probability through U(x) or the time delay associated with the 
same potential) are registered by receiver. Since the scattering properties of 
the fluxons depend on U(x) that is controlled by the measured system, the 
scattered fluxons contain the information about the state of the system. The 
injection frequency is set sufficiently low so that only one fluxon at a time 
moves inside JTL. 

surement is to detect the variations in the junction impedance caused by the measured 
system. Such " impedance -measurement" schemes j^J enable one to perform measure- 
ments continuously, but they typically require narrow-band coupling between the detector 
and the measured system which in turn limits their time resolution. 
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Figure 6.1: (a) Equivalent circuit of the flux detector based on the Josephson 
junction transmission line (JTL) and (b) diagram of scattering of the fluxon 
by the potential U(x) that is controlled by the measured qubit. The fluxons 
are periodically injected into JTL by the generator and their scattering char- 
acteristics (transmission and reflection coefficients t(k), r(k)) are registered by 
the receiver. 



In what follows, we are interested in the regime of low segment inductances 
L <C &o/I c where $o = Tth/e is the magnetic flux quantum. This-way opti- 
mized JTL acts as a uniform ballistic fluxon channel. The phase difference 
across each segment is small, and the junction can be viewed as distributed 
structure described by standard sine-Gordon Lagrangian, 

L = ^JT L j dx ji(<9 r 0) 2 - i(a^) 2 + cos(0) - <P (e) (x)d x <P^ , (6.1) 

where x is expressed in the units of Josephson penetration depth A j = ^le^jej. 
The Josephson energy is replaced by Josephson energy density ej = $ J c /2tt 
(i.e. I c — > J c ), while the segment inductance L becomes junction inductance 
density / which gives inductive energy = ($ /27r) 2 /2Z units of "energy" x 
"distance". The time is conveniently expressed in the units of inverse plasma 
frequency cj 2 = €j/c where c is the capacitance density of the line. The 
commutation relation for the charges and phases of individual junction give 
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the following equal-time commutation relations for the field 0(r, x): 

[<J ) (x),d T( j ) (x')}=p 2 5(x-x') J (6.2) 

where (3 2 = (4e 2 / K) / measures the wave impedance \/Iq/c of the JTL 
in the absence of Josephson tunneling relative to quantum resistance. 

Known results for quantum sine-Gordon model (e.g. see jHHj) show that 
for (3 2 > 87r, i.e. a//o/ c o > h/e 2 ~ 25/cf2, quantum fluctuations of the field <fi 
completely destroy the quasiclassical excitations of the junction. This situation 
is similar to that in small Josephson junction jM]; and the dynamics of the 
supercurrent flow in JTL should then be described in terms of tunneling of 
individual Cooper pairs [T3]. While this regime might be reachable in very 
narrow JTLs of sub-micron width [HZ] we assume in this work the more typical 
situation when the impedance v^o/ c o is on the order of 100f2, and (3 2 <C 1. In 
this case, the JTL supports a number of quasiclassical excitations including 
most importantly for this work topological solitons that each carry precisely 
one quantum of magnetic flux. The dynamics of such "fluxons" is equivalent 
to that of stable, in general relativistic particles [33] with terminal velocity 
c = XjLUp and mass 

m = Vl^/XWp- (6-3) 

Another type of quasiclassical excitations in JTL are the small-amplitude plas- 
mon waves with frequency 

u{k) = {u 2 p + c 2 k 2 f/ 2 , (6.4) 

where k is the plasmon wave vector. 

In this work we are interested in "non- relativistic" regime of fluxon dynam- 
ics, when the fluxon velocity u is small in comparison to terminal velocity c, 
i.e. nCc. The equations (|6.3|) and (J6.4j) show that in this regime, the fluxon 
terminal energy e = k 2 /2m can be made smaller than the lowest plasmon 
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energy hut. 



> P - 



(2u/cP) 2 , 



(6.5) 



so that for u < c(3 the fluxon cannot emit a plasmon even when it is scattered 
by non- uniformities of the JTL potential |98| . Intrinsic dissipation associated 
with emission of plasmons is then suppressed, and the fluxon motion in JTL is 
elastic, provided that other, "extrinsic", sources of dissipation are sufficiently 
weak. Another case when plasmon dissipation is reduced including the cases of 
fast fluxons, is a situation of sufficiently smooth JTL potential. In this case the 
plasmon emission precesses are suppressed even when they are energetically 
possible. Although the JTL operation as the flux detector should be possible 
for any strength of fluxon dissipation, significant dissipation would prevent the 
detector from reaching the quantum-limited regime. 

The shape of the scattering potential U (x) for fluxons created by the mea- 
sured system is determined by the convolution of the distribution of the flux 
$( e ) with the distribution of current in each fluxon P^lllUUj . an d can be written 
as: 



where 4>(x) is the shape of the fluxon that in general can be distorted by the 
potential U(x) itself. If the potential is smaller than u p , or does not vary 
appreciably on the scale on the size of the fluxon given by Aj, the changes in 
the fluxon shape are negligible and one can use in equation ()6.6|) the regular 
fluxon shape in the uniform case. In non-relativistic limit the shape of the 
fluxon is given by sine-Gordon kink soliton 4>o(x) = 4 arctan[exp(x/Aj)]. One 
of the implications of equation ()6.6|) for the JTL detector is that the width of 
the scattering potential U(x) cannot be made smaller than Aj. 

If the measured system is not coupled to JTL inductively as in figure 16.11 
but rather galvanically, and injects the current with density j^ e '(x) in the nodes 
of JTL as show in in figure 16.21 the potential created in the junction is still 




(6.6) 
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Figure 6.2: Equivalence between the magnetic and galvanic coupling to the 
detector. External fluxes $1^ through the JTL cells induce the circulating 
currents j n in them, which in turn create the currents jn = j n ~ jn-i through 
the JTL junctions. The situation would be the same if the currents jn^ are 
injected directly into the JTL junctions by external system. 



given by the equation ()6.6|) if one substitutes 

(l// )<9$ (e) /<9x = j ie \x). (6.7) 

Another important remark is that although the assumed condition (3 2 <^ 1 
makes quantum fluctuations of the fluxon shape small, the dynamics of the 
fluxon as whole can still be completely quantum, as recently observed experi- 
mentally |lUlj . Achieving the limit of small (3 2 demands low capacitance of the 
junction. This can be achieved by in reality constructing the JTL consisting 
of small junctions and inductances, rather than out of uniform electrodes (see 
figure l6~5j) . As discussed in |lU2j . the above description is still valid if the cell 
size is made smaller than the fluxon length Aj. 

Operation of the JTL detector is based on free, non- relativist ic dynamics 

of fluxons of mass ()6.3|) on the potential described by ()6.6|) and (|6.7|) . It also 

requires that the fluxons are injected by the generator into one end of the 

JTL and detected by the detector on the other end ( figure l6~Tj) . Both of these 

circuits can be designed following the general principles of SFQ electronics 

|lU3j and will not be explicitly discussed here 2 . It is simply assumed that the 

ends of JTL are matched appropriately to these circuits so that the fluxons 

2 Adaption of SFQ electronics to qubit applications [1041 1105| is an important and not 
yet fully understood part of the development of scalable superconducting qubits 
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Ground plane 



Figure 6.3: Geometry of one cell of the multi-layer JTL: top view and verti- 
cal cross-section. Internal layers of the structure are used to create Joseph- 
son junctions (JJs), while the two external layers, the ground plane and the 
counter-electrode, separated vertically by the distance U, define inductance L. 
The length a of the cell is set by the distance between the junctions. 

can enter or leave the JTL without reflection and are injected into the JTL in 
appropriate quantum state. 

Initial quantum state ip(x,t = 0) of an injected fluxon is characterized by 
the average fluxon velocity Uq and the shape of the wave packet ipa(x) defining 
its position: 

if,(x, t = 0) = ip (x)e ikoX , h = tjiuq. (6.8) 

As it will be seen form the discussion in the next section, many properties of the 
JTL detector are independent of the specific shape of the wave packet ip ( x ), 
as long as it is well localized both in coordinate and momentum space of the 
fluxon. These properties depend on the wave packet width w in coordinate 
space and the corresponding width 8k ~ 1/w in momentum space. These 
parameters should satisfy two obvious conditions: (1) w <C I, where I is the 
total length of JTL, and (2) 8k <C k$. We assume a stronger version of the 
latter condition that follows from the requirement that the broadening of the 
wave packet by 8x ~ 8kt/m due to momentum uncertainty and travel time 
t ~ 1/uq is negligible in comparison to the original width: 8x <C w. The 
discussion below shows that this requirement is necessary for the quantum- 
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limited operation of the detector in the time-delay mode. The two conditions 
put together limit the width w such that: 

(l/k ) 1/2 < wl. (6.9) 

In the cases when it is necessary to specify the shape of the wave packet 
ipo( x ) we wm t & k e it to be Gaussian: 

i^q{x) = {Tiw 2 )~ l l A exp 

where x is the initial fluxon position in the JTL. Besides being well local- 
ized both in coordinate and momentum space, the wave pocket (j6.10j) can be 
obtained as the result of fluxon generation process naturally implementable 
with SFQ circuits. The fluxons satisfying (J6.9)) and (jfi.lOj) can be obtained in 
two steps: (1) relaxation of the fluxon to the ground state of weakly damped 
and nearly quadratic Josephson potential with required width w of the wave 
function of the ground state, and (2) rapid acceleration to velocity u after the 
potential is switched off. The fluxons injected in JTL in the state (|6.8|) can be 
used for measurements as described in the following section. 

6.2 Measurement Dynamics of the JTL De- 
tector 

The measurement by the JTL detector consists qualitatively in scattering of 
the fluxons in the JTL of the potential controlled by the measured system. 
This process has the simplest dynamics if the measured system is stationary. 
In this case, it is convenient to consider it in the basis of the eigenstates \qj) of 
the system operator (e.g. magnetic flux in the qubit loop in the example shown 
in figure I6~T]) which couples to the JTL. In each state \qj) the system creates 
different potential Uj(x) for the fluxons propagating through JTL. Different 
realizations of Uj(x) produce different scattering coefficients for injected flux- 



(x — x) 
2w 2 



2 



(6.10) 
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ons: the amplitude tj(k) of transmission and the amplitude rj(k) of reflection 
back to toward the generator. Both amplitudes depend on fluxon momentum 
k. Since the scattering amplitudes depend on the state \qj) of the measured 
system, scattered fluxons carry information about \qj). 

For the JTL detector, the detector-system entanglement arises as a result 
of fluxon scattering, and the rates of information acquisitions and back-action 
dephasing can be expressed in terms of the scattering parameters. In this 
respect JTL is very similar to QPC detector and its characteristics can be 
obtained following the similar derivation |lU2j for QPC. 

Evolution of the density matrix p of the measured system in scattering of 
one fluxon can be obtained by considering first the time dependence of the 
total wave function of the fluxon injected into JTL and the wave function 
J2j a jl<?j) °f the measured system: 



Here the initial fluxon wave function ip(x,t = 0) is given by equation (|6.8|) and 
its time evolution ip(x,t) depends on the realization Uj(x) of the scattering 
potential created by the measured system. Evolution of p is then obtained by 
taking the trace over the fluxon part of the wave functions ([6.11)1 : 



Qualitatively, the time evolution in equation ([6.11)) describes the propaga- 
tion of the initial wave packet (|6.8|) towards the scattering potential and then 
separation of this wave packet in coordinate space into the transmitted and 
reflected parts that are well-localized on the opposite sides of the scattering 
region. If we assume that the scattering potential U (x) has a simple shape 
(e.g. without narrow quasi-bound states) and is non-vanishing only in some 
small region of size a <C /, the time t sc from the fluxon injection to completion 
of the scattering is not much different from the time 1/uq of free fluxon prop- 



(6.11) 



j j 




(6.12) 
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agation through JTL. Then at time t > t sc , the separated wave packets move 
in the region free from j'-dependent scattering potential and the unitarity of 
the quantum mechanical evolution of ipj (x, t) implies that the overlap of the 
fluxon wave functions in (J6.12)) becomes independent of t. 

This overlap can directly be found in momentum representation: 



dxi>i(x,t)^*(x,t) = / dk\b{k)\ 2 [Ut*+r ir *] . 



(6.13) 



Here b(k) is the probability amplitude for the fluxon to have momentum k in 
the initial state (|6.8[) . e.g. in the case of the Gaussian wave packet ()6.10|) 



b(k) = {w 2 /n)- 1/A 



exp 



(k - k fw 2 



i{k — ko)x 



(6.14) 



Equations (J6.12)) and (J6.13j) show that the diagonal elements of the density 
matrix p do not change in the process of scattering of one fluxon, while the 
off-diagonal elements are suppressed by the factor 



V 



dk\b(k)\ 2 [uWtjWih) + r i (k)r* j (k)] 



< 1 



(6.15) 



where the inequality is the consequence of Schwartz inequality as mentioned 
in the previous chapter. After summing the suppression factors for all fluxons 
injected at a frequency /, the rate of back-action dephasing is obtained as 



-/In 



dk\b{k)\ 2 \ti{k)tj{ky + n(k)r*(k)] 



(6.16) 



Equation (|6.16|) is similar, but not identical, to the back-action dephasing 
rate by the QPC detector. The main difference is at the stage of momen- 
tum summation, since different momentum electrons scatter independently, 
while this is not the case for fluxons. They are constrained by the condition 
that whole fluxon is either transmitted or reflected by the potential, since the 
scattering events of different fluxons are well separated in time. As we will 
see below, this difference makes it possible to operate the JTL detector in 
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time-delay mode that is not possible with the dc-biased QPC detector. 

As examples of applications of equation (|6.16J) we consider several specific 
cases motivated by the measurement regimes discussed below. In one, the 
phases of the scattering amplitudes are assumed to cancel out form equation 
(I6.16JI . while the variation of the absolute values of the amplitudes is small. 
Then, the dephasing rate (j6.16|) can be expressed in the terms of variations 
5Tj(k) of the fluxon transmission probability in different states \j) around 
some average transmission T(k): \tj{k)\ 2 = T{k) +STj(k), STj(k) <C T{k). To 
the lowest non- vanishing order in STj(k) we get: 



This equation describes "linear response" regime of operation of the JTL de- 
tector, when its properties follow from the general theory of linear detectors 
|12| . In particular the dephasing (J6.17|) can be understood as being caused 
by the back-action noise arising from the randomness of the fluxon transmis- 
sion/reflecion process. 

In the "tunnel" limit of weak transmission |ij(A;)| <C 1, and still under the 
assumption of cancelling scattering phases, the equation (jfi.lfij) reduces to: 



The last two equations have corresponding analogues in the case of QPC de- 
tectors una una El- 

As the last example, that does not have an analogue in the QPC physics, 
we consider the situation when the reflection is negligible, |rfc(t)| = 0, and 
the system- JTL interaction modifies only the phases Xj(fc) of the transmission 
amplitudes: tj(k) = e l% ^ k \ so that 



In coordinate representation this means that the scattering only affects the 




(6.17) 




(6.18) 




(6.19) 
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position and the shape of the transmitted wave packet. If the width 5k of 
initial fluxon state in the momentum representation is sufficiently narrow, and 
the phases Xj{k) do not vary strongly over the momentum range, they can be 
approximated as 

Xj(k) = Xj(ko) -(k- ko)x,, Xj = x'j(k )- (6.20) 

This approximation neglects the distortion of the wave packet's shape during 
the scattering process, while it does take into account the potential induced 
shift Xj along the coordinate axis. This shift can directly be related to the 
"time delay" Tj = Xj/uq due to scattering. The value of the delay can take 
on positive and negative values depending on the scattering potential. The 
dephasing rate equation ()6.19|) is conveniently expressed in coordinate repre- 
sentation in the terms of initial fluxon wave packet ipo(x): 



T i:j = - fin 



dxi/)o(x — Xi)ipl{x — Xj) 



(6.21) 



This expression explicitly shows that the back-action dephasing by the JTL 
detector arises form the entanglement between the measured system and the 
scattered fluxons which are shifted in time by an interval dependent on the 
state of the measured system. The degree of suppression of coherence between 
the different states of the measured system in scattering of one fluxon is deter- 
mined then by the magnitude of the relative shift of the fluxon in these states 
on the scale of wave packet width. For instance, if the wave packet of initial 
fluxon is Gaussian (|6.1U|) the equation (|6.21|) gives: 

r« = f i - E ^r 1 - (6-22) 

Back-action dephasing represents only a part of the measurement process. 
The other part is the information acquisition by the detector about the state 
of the measured system. In the case of the JTL detector, this information is 
contained in the scattering characteristic of fluxons, and the rate of its acqui- 
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sitions depends on the specific characteristics recorded by the fluxon receiver. 
There are at least two possibilities in this respect. One is to detect the prob- 
ability of fluxon transmission through the scattering region (or, equivalently, 
the corresponding probability of the fluxon reflection back to the generator). 
Another possible detection scheme is for the receiver to measure the time delay 
associated with the fluxon propagation through the JTL. Even if the measured 
system changes the potential Uj(x) in such way that the fluxon transmission 
probability is not affected, the potential can still change the fluxon's propaga- 
tion time, which will contain the information about the state of the measured 
system. In general, one can have a situation when the information is con- 
tained both in the changes of propagation time and transmission probability, 
and in order not to lose any information one would need to detect both scat- 
tering characteristics. In this work, we consider only the two "pure" cases of 
transmission and time-delay detection modes. 

6.2.1 Transmission Detection Mode 

If the fluxon receiver of the JTL detector records only the fact of fluxon arrival 
at the recever, then only the modulation of the fluxon transmission probability 
by the measured system conveys information about the system. The informa- 
tion contained in all other features of the scattering amplitudes (e.g. their 
phases [TJ] or propagation time) is lost in the receiver. In this case the rate of 
information acquisition can be calculated simply by starting with the proba- 
bilities of transmission/reflection Tj and Rj, when the measured system is in 
the state \qj): 



The independence among the different fluxon events implies that the prob- 
ability to have n transmitted out of N incident fluxons is given by binomial 
distribution 




(6.23) 



Pj (n) = CftTj'Rj 



(6.24) 
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The task of distinguishing different states \qj) of the measured system is trans- 
formed into distinguishing the probability distribution Pj(n) for different j's. 
If the ffuxons are injected into the JTL periodically with frequency /, the 
number of the scattering events increases with time, N = ft, and the proba- 
bility distribution f)6.24j) becomes more and more strongly peaked around the 
average number NTj of transmitted fluxons. Thus the states with different 
Tj are distinguishable with increasing certainty. The rate of the increase of 
this certainty can be characterized quantitatively by some overlap of different 
probability distributions Pj(n). While there are different ways to characterize 
the overlap of different probability distributions [T^], the characteristic which 
is appropriate in quantum measurement context [HJ |SE] is closely related to 
"fidelity" in quantum information [T3]: ^ n [pj(^)Pj(^)] 1//2 - The rate of infor- 
mation acquisition is then naturally defined as 

Wn = -l/thxY,\Pi{n.)Pj{n)] 1/2 , ( 6 - 25 ) 

n 

which after substitution of (J6.24J1 becomes: 



Wij = -fhx 



y/Wj + y/R^\, (6-26) 



where the transmission and reflection probabilities are given in (|6.23j) . 

The equation (|fi.2fij) characterizes the information acquisition rate of the 
JTL detector in transmission detection mode. For an arbitrary detector the in- 
formation acquisition is smaller or equal to the back-action dephasing, and the 
regime when the two rates are equal is quantum-limited. Comparing equations 
(Ifi.lfij) and (J(L26|) one can see that for our detector, indeed 

Wij < Iy, (6.27) 

and that the equality golds if several conditions are satisfied. First two con- 
ditions require that there is no information in the phases of the transmission 
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amplitudes: 

<j>i(k) = <t>j{k), (j)j(k) = argitjW/rjik)) (6.28) 
Xj(k) — Xi{k) = const. (6.29) 

Equation ()6.29)) means that the difference of the phases Xj{k) of the transmis- 
sion amplitudes should be effectively independent of the momentum k in the 
relevant range set by the fluxon distribution |fr(A;)| 2 over k. The two conditions 
()6.28|) and ()6.29|) have different physical origin. The first implies that the scat- 
tered states contain no information on \qj) that can be obtained in principle by 
arranging an interference between the transmitted and reflected parts of the 
wave function ^2] • in the practical terms, the safest way to satisfy this require- 
ment is to make the scattering potential symmetric, i.e. Uj(—x) = Uj(x) for 
all states \qj). The unitarity of the scattering matrix for the fluxon scattering 
in the JTL implies in this case that <pj = ir/2 for all j's. 

Condition ()6.29|) means that no information on \qj) is contained in the 
shape or position of the transmitted wave packet that would be lost in the 
fluxon receiver operating on the transmission-detection mode. The similar 
condition for the reflected fluxon is implied by both requirements ()6.28|) and 
()6.29|) . In general, condition (J6.29J) requires that the spread of the initial wave 
packet Sk in the momentum space gives the rise to the energy uncertainty of 
the fluxon 5e ~ uo5k, that is much smaller than the energy scale Q of the 
transparency variation of the scattering potential Uj(x). 

One more condition of the quantum-limited operation of the JTL detector 
in the transmission-detection mode is that the fluxon transmission probabilities 
are effectively momentum and therefore energy independent in the respective 
relevant ranges: 

Hk)\ 2 =T r (6.30) 

This condition requires that 5e <C Q. It is more restrictive than the corre- 
sponding condition for the QPC detector which can be quantum-limited even 
in the case of energy-dependent transmission probability O I1U8| [L09J. To 
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obtain equation (J6.30J) one starts from the back-action dephasing rate (j6.16|) 
which can be written as 

Iy = -fin J dk\b(k)\ 2 [\t t (k)t*(k)\ + \ ri (k)r*(k)\] (6.31) 

under the conditions ()6.28|) and ()6.29|) . Schwartz inequality for the functions 
Uk)\: 

^[W 3 > J dk\b(k)\ 2 \U(k)tj(k)l 

and the similar inequality for |rj(/c)| show that from ()6.3H) and the in- 
formation acquisition rate Wij satisfy the inequality (|6.27|) . Furthermore, the 
equality is reached only when 

\t j (k)\=£ j t(k), \r j (k)\=? j r j (k), (6.32) 

and the ratios \U{k)\/\tj{k) \ and \ri(k)\/\rj(k)\ are independent of k. Similarly 
to (|6.29|) equations (|6.32j) demand that no information about the state of the 
measured system is contained in the shape of transmitted or reflected wave 
packets. In general, when both transmission and reflection probabilities are 
not small the two relations in ()6.32|) are incompatible. They have only a 
trivial solution in which all amplitudes are independent of k in the relevant 
range of fiuxon momentum distribution, thus proving the equation ()6.30|) for 
Tj ~ Rj ~ 1/2. The transmission and reflection probabilities have the roughly 
the same magnitude when the fiuxon energy e is close to the maximum U of 
the scattering potential. In this case, small spread of e: 5e <C Q implies that 
the characteristic range of the scattering potential should be small: a <C w. 

The condition (|6.30J) of the quantum-limited operation is not necessary 
when either Tj < 1 or fij < 1. In this case one of the relations in ()6.32|) 
applies to the linear response regime, when the variations 5Tj of JTL trans- 
parencies between the different states \qj) are small and the back-action de- 
phasing is given by (J6.17|) . Expanding (|6.26j) in 8Tj : Tj = T + 5Tj , where all 
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transparencies are defined as in (|6.23|L we get: 



IT(1 — T) 



W v = f 0Z1 it ■ (6-33) 



This equation differs from the equation (|6.17|) only by the order in which the 
integration over the momentum is performed. This means that the information 
acquisition and dephasing rates satisfy in general the inequality (J6.27j) and are 
equal only if the transparency of JTL is constant in the relevant momentum 
range. 

In the linear response regime, each individual ffuxon carries only small 
amount of measurement information, and it is convenient to employ the quasi- 
continuous description in which the fluxon receiver acts as the voltmeter reg- 
istering not the individual fluxons, but the rate of arrival of many fluxons 
represented by voltage V(t) across the segment junctions of JTL. The average 
voltage for state \qj) is 

(V(t)) = fT&o, (6.34) 

where (••■) implies average over the scattering outcomes and time t of the 
fluxon injection cycle. Because of the randomness of the fluxon scattering, the 
actual voltage fluctuates around the average value (J6.34J) even at low temper- 
atures. These fluctuations can be attributed to the shot noise of fluxon and 
its spectral density 

S v (uj) = [ dre'^ ((V(t + r)V(t)) - (V(t)} 2 ) , (6.35) 



which is constant at frequencies uo below the fluxon injection frequency /. 
Straightforward calculation similar to that for regular shot noise shows that 
the constant value of the spectral density is 

S = fT(l-T)& , (6.36) 

where in the linear response regime we can neglect the small differences STj 
of transparencies between the different states \qj) in the expression for the 
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noise. The equation shows that in accordance with the general theory of linear 
quantum measurements |12j . the information rate ()6.33|) as the rate with which 
one can distinguish dc voltage values ()6.34jl in the presence of white noise with 
spectral density (|6.36|) . 

To conclude this subsection, we arrive the the same conclusions by following 
the conditional evolution approach. This is easily accomplished by looking at 
the evolution of the density matrix elements in the conditional evolution as 
given in (J5.17|) and (J5.18j) . After applying the requirement that the qubit stays 
"pure" during the ideal measurement (i.e. pooPn = |/°oi| 2 ) the cross-integrals 
of the transmission/reflection rates need to satisfy: 



There two conditions are only satisfied if the rates tj,rj are momentum inde- 
pendent in the range of the injected fluxons and if the phases of the rates are 
the same. This is equivalent to the conditions (|6.28|) and (J6.32J) . which leads 
to the main conclusion of this section that in transmission-detection mode 
the most relevant regime Tj ~ 1/2 of maximal detector response to the input 
signal, the conditions of quantum-limited operation are given by conditions 
ijOSt and IjOty . 

6.2.2 Time-Delay Detection Mode 

Since the range a of the scattering potential ()6.6|) cannot be smaller than 
the fluxon size Aj it might be difficult in practice to satisfy the a <C w con- 
dition needed for the quantum-limited operation of the JTL detector in the 
transmission-detection mode. For quasiclassical potential barriers that are 
smooth on the scale of the fluxon wave packet, a > w, the "transition" region 
of energies near the top of the barrier, where the reflection and transmission 
amplitudes have comparable magnitude, is narrow. If the interval 8e of the 




1. 



1. 
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energies avoids this narrow region, then either the transmission or reflection 
coefficient can be neglected. Ballistic motion of the fluxons in this region 
still contains information about the potential Uj(x) that can be used for mea- 
surement. This information is contained in the time shift Tj caused by the 
propagation though the region with non-vanishing potential rather than in 
transmission/reflection probabilities. Quantum mechanically, the time-shift 
information is contained in the phases of the scattering amplitudes. To be 
specific, we discuss here the regime where the JTL detector is operated in 
time-delay mode using the transmitted fluxons, = 1- I n the energy 

range where !?".,■(&) | = 1, the same detection process is possible using the re- 
flected fluxons. The only advantage of the transmission case is the possibility 
to use the fill range of values of the scattering potential: Uj(x) ^ 0. 

For sufficiently smooth Uj(x) the phase Xjify of the transmission amplitude 
can be calculated in quasi-classical approximation: 



If the potential is weak, Uj(x) <C e, the potential-induced contribution to the 
phase (|6.37j) can be expressed as 



Under the adopted assumption of quasiclassical potential and vanishing 
reflection, the condition ()6.9j) of the negligible broadening of the fluxon wave 
packet during free propagation through the JTL is still walid in the presence of 
the potential. In this case one can use the approximation (|6.20|) for the phases 
Xj(k) which implies that the wave packet shifts as whole and without distor- 
tion. The potential induced part of the shift Xj = — x'j{ko) can be obtained 
from the equation (|6.37|) and has the classical form: 




(6.37) 




(6.38) 




(6.39) 
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For weak potential, Uj(x) <C e: 



Xj = / dxUAx 

3 2eJ jK 



(6.40) 



Back-action dephasing rate by the JTL detector in this regime is given by 
equations ()6.21)1 and ()6.22|) . The information about the states \qj) contained 
in the shift Tj of the fluxon in time or equivalently in coordinate Xj = uoTj, can 
be read-off by distinguishing different shifts Xj against the background of finite 
width w of the fluxon wave packet ipo(x). Since |^o(20| 2 gives the probability 
of finding the fluxon at coordinate x, this task is equivalent to the task of 
distinguishing two shifted probability distributions (see figure I6.4|) that was 
discussed above for the transmission-detection mode. Similarly to equation 
(I6.25j) we can write the information acquisition rate of the JTL detector in the 
time-delay mode as follows: 



Comparing this to the dephasing rate (j6.21|) . we see that in general the two 
rates satisfy the inequality (J6.27J) as they should. The rates are equal if the 
phases of the initial packet ipo(x) of the injected fluxon is independent of x, 
i.e. if ipo(x) is real. In particular case of the Gaussian wave packet, the 
detector is quantum-limited (Wij = r^), and the two rates are given by ()6.22|) . 
These considerations also imply that the JTL detector in the time-delay mode 
would loose the property of being quantum-limited if the fluxon wave packet 
spreads noticeably during the propagation through JTL. This process creates 
non-trivial x-dependent phase of the wave packet and makes the information 
acquisition rate smaller than the back-action rate. 

The derivation of the conditional evolution equations introduced in the last 
chapter can be repeated with only minor modifications in time-delay mode of 
detector operation. In this case, the different classical outcomes of measure- 
ments are the observed instances of time when the fluxon reaches the receiver, 
that for convenience can be directly translated into different fluxon positions 




(6.41) 
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Figure 6.4: Illustration of the information acquisition process by the JTL 
detector in the time-delay mode. The fluxon wave packet is shifted by the 
distance Xj dependent on the state \qj) of the measured system. In conditional 
description, observation of the fluxon with position x' changes the system wave 
function according to equation (|6.42l) . 

x at some fixed time. If the fluxon is observed at the position x' in a given 
injection cycle ( figure l6~4~j) . the evolution of the amplitude oij of the measured 
system due to interaction with the fluxon is: 



Qualitatively, and similarly to the conditional evolution in the transmission- 
detection mode, the sequence of transformations f)6.42j) describes weak mea- 
surement: gradual localizationof the system wave function to one of the states 
\j) with increasing number of fluxon scattering events which lead to the grad- 
ual accumulation of the information about the system. In contrast to the 
transmission-detection mode, conditional evolution mode (I6.42J1 always re- 
mains pure under the symplifying assumptions adopted in this work for the 
JTL detector: coherent propagation of fluxons with wave packet of fixed form 
i/jq(x) injected into the JTL with the same energy. 




(6.42) 
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6.3 Non-Quantum-Limited Detection 



Quantum-limited operation of the JTL detector discussed in the preceding 
sections requires quantum coherent dynamics of fluxons in the JTL. While this 
dynamics can be observed experimentally |lUlj . the task of realizing it is very 
difficult. From this prospective it is important that several attractive features 
of the JTL detector, e.g. large operating frequency and reduced parasitic 
dephasing during the time intervals between the fluxon scattering, remain even 
in the non-quantum-limited regime. Our discussion above was limited on the 
fundamental back-action dephasing which in its quantum limit is unavoidable 
part of measurement process. Although the Josephson junctions of JTL and 
those of "external" parts of the JTL detector (generator and receiver) can 
give rise to dissipation and dephasing not related directly to measurement, in 
the JTL geometry (figure I6.1J) parasitic dissipation is suppressed due to the 
screening by the supercurrent flow in the JTL junctions |1U4[ 111)5] . 

The dominant deviation from the quantum limited detection should be as- 
sociated then with the fluctuations in the fluxon motion. These fluctuation 
make the dephasing factor r] (|6.15|) due to scattering larger than the amount of 
the information conveyed by the fluxon scattering. Some interesting measure- 
ment strategies are still possible with the JTL non-ideality of this type. The 
most natural example is QND measurement briefly discussed in the previous 
chapter (e.g. see O EH]) which in principle make the detector back-action 
irrelevant. The possibility of timing the fluxon scattering in the JTL detector 
makes it particularly suitable for the "kicked" version of the QND qubit 
measurement. 

Consider the qubit Hamiltonian 

H = ~^a x (6.43) 

which performs quantum coherent oscillations with frequency A . The qubit 
is coupled to the JTL through its a z operator, i.e. the states \qj) used in the 
discussion in the preceding sections, are the two eigenstates of a z . We assume 
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Figure 6.5: Schematics of the QND fluxon measurement of a qubit which 
suppresses the effect of back-action dephasing on its coherent quantum oscil- 
lations. The fluxon injection frequency / is matched to the qubit oscillation 
frequency A: / ~ A, so that the individual acts of measurement are done 
when the qubit density matrix is nearly diagonal in the a z basis, and the 
measurement back-action does not introduce dephasing in the oscillation dy- 
namics. 

that the Hamiltonian (J6.43)) already includes renormalization of parameters 
due to the qubit-detector coupling. If the qubit oscillations are weakly de- 
phased at the rate 7 < A (e.g., by residual parasitic dissipation in the JTL 
detector), the time evolution of the qubit density matrix p during the time 
intervals between the successive fluxon scattering processes can be written in 
the following form in the <r 2 -representation: 



where r(t = 0) = ±1 depending on weather the qubit starts at t — from the 
o z = 1 or a z = — 1 state. 

The fluxon scattering at times t n = n/f leads to suppression of the off- 




(6.44) 



f = — iAr, r = x + iy, 



(6.45) 
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diagonal elements of p: 



y{t$) =vy(t n )- 



(6.46) 



However, if the ffuxons are injected in the JTL with the time interval 1/ / close 
to the half-period ir/A of the qubit oscillations f figure l6~o"|) . the qubit density 
matrix (|6.45J) is nearly diagonal at the moments of fluxon scattering, y(t n ) <C 1, 
and suppression ()6.46J) does not affect strongly the qubit oscillations. 

Such a QND qubit measurement is possible with the JTL detector op- 
erating in either detection mode; to be specific we assume the transmission 
mode. The dependence of the fluxon transmission probability on the qubit 
state can be written then as T+a z 5T . We consider the linear-response regime, 
5T <C T, when the scattering of one fluxon provides only small amount of in- 
formation about the qubit oscillations. For quantum-limited detection, the 
linear-response condition 5T <C T implies that r] — > 1. In the non-quantum- 
limited case, the back-action can be stronger, and we take i] to be an arbitrary 
factor within the [0, 1] interval. 

If the frequency / is matched precisely to the qubit oscillations, / = A/tt, 
the matrix ()6.45|) is diagonal at the times of measurement, y(t n ) = 0, and 
detector does not affect at all the qubit dynamics. If the mismatch is non- 
vanishing but small, 5 = A/ f — 7r <C 1, diagonal elements of p (J6.45)) evolve 
quasi- continuously even if suppression factor 7] is not close to 1. Equations 
()6.45|) and ()6.46|) give the following equation for this quasi-continuous evolu- 
tion: 



In the assumed linear-response regime, the qubit oscillations manifest them- 
selves as a peak in the spectral density Sy(u) (|6.35|) of the voltage V across 
the JTL junctions. For / ~ A/ir the oscillation peak in Sv{uj) is at zero fre- 
quency. Equations (|6.44|) and (|6.47J) describing the decay of correlations in the 
qubit dynamics in the a z basis imply that the oscillation peak has Lorentzian 




(6.47) 
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shape: 



b v {uj) = bo^ 2 , r2 — , (6.48) 

and the oscillation line width T is 

For the quantum-limited detection, r\ — > 1, the equation ()6.49|) reproduces 
previous results for the QND measurement [13] if one introduces the back- 
action dephasing rate 1^ = f(l — rj). In this case, equation ()6.49|) is valid for 
sufficiently small mismatch between the measurement and oscillation frequen- 
cies, \8\ <C 1— 77. For larger S, the oscillation peak in the detector output Sy(u)) 
moves to finite frequency A — irf and the QND nature of the measurement 
is lost [H]. Equation ()6.49|) shows also that in the limit of "projective" mea- 
surements f] = 0, broadening of the oscillation peak is weaker, and the peak 
remains at zero frequency for all reasonable values of the detuning parameter 
1 5 1 1. Therefore, the stronger back-action of the JTL detector is advan- 
tageous for the QND measurements of coherent qubit oscillations. The last 
remark is that although our discussion here assumed that the fluxon arrival 
times are spaced exactly by 1/ /, equations (|6.48|) and (|6.49|) should remain 
valid even in the presence of small fluctuations of the measurement times. 
These fluctuation can be described by taking into account that the detuning 
\S\ cannot be made smaller that the relative line width of the fluxon generator. 

This concludes the discussion about the ballistic JTL detector for the flux 
qubits - the last part of this thesis. 



106 



Bibliography 



[1] Allan Turing. Proc.Lon.Math.Soc, 42:230, 1936. 
[2] A. Church. Am. J. Math., 58:345, 1936. 

[3] R. Landauer. Irreversibility and heat generation in the computing pro- 
cess. IBM J. of Res. and Devel, 3:183, 1961. 

[4] C. Bennett. Logical reversibility of computation. IBM J. of Res. and 
Devel, 17:525, 1973. 

[5] D. Deutsch. Quantum theory, the church turing principle and the uni- 
versal quantum computer. Proc. Roy. Soc. London Ser. A, 400:73, 1985. 

[6] D. Deutsch and R. Jozsa. Rapid solution of problems by quantum com- 
putation. Proc. Roy. Soc. London Ser. A, 439:553, 1992. 

[7] D. Simon. On the power of quantum computation. In Proceedings 35th 
Annual Symposium on Foundations on Computer Science, Los Alamitos, 
CA, 1994. IEEE Press. 

[8] P. Shor. Algorithms for quantum computation: discrete logarithms and 
factoring. In Proceedings 35th Annual Symposium on Foundations on 
Computer Science, Los Alamitos, CA, 1994. IEEE Press. 

[9] L.K. Groover. Quantum mechanics helps in searching for a needle in a 
haystack. Phys. Rev. Lett., 79:325, 1997. Also in quant-ph/9706033. 

[10] D. DiVincenzo and B. Terhal. Phys. World, 11(3):53, 1998. 



107 



[11] H.B. Callen and T.A. Welton. Phys.Rev., 83:34, 1951. 

[12] D.V. Averin. Quantum Noise in Mesoscopic Physics, volume 97, page 
229. Kluwer Academic Publishers, 2003. also cond-mat 0301524. 

[13] Michael A. Nielsen and Isac L. Chuang. Quantum Computation and 
Quantum Information. Cambridge University Press, 2000. 

[14] D.V. Averin, A.B. Zorin, and K.K. Likharev. Zh.Exp.Theor.Fiz., 88:692, 
85. (Also Sov.Phys. JETP 61,407). 

[15] M. Biittiker. Phys.Rev. B, 36:3548, 1987. 

[16] D.L. Flees, S. Han, and J.L. Lukens. Phys.Rev. Lett, 78:4817, 1997. 

[17] V. Bouichiat, D. Vion, P. Joyez, D. Esteve, and M. Devoret. J. of Super- 
cond., 12:789, 1999. 

[18] Y. Nakamura, Yu.A. Pashkin, and J.S. Tsai. Nature, 398:786, 1999. 

[19] D. Vion, A. Assime, A. Cotter, P. Joyez, H. Pothier, C. Urbina, D. Es- 
tave, and M. H. Devoret. Science, 296:886, 2002. 

[20] Yu. A. Pashkin, T. Yamamoto, O. Astafiev, Y. Nakamura, D.V. Averin, 
and J.S. Tsai. Nature, 421:823, 2003. 

[21] A. Wallraff, D. I. Schuster, A. Blais, L. Frunzio, R.-S. Huang, J. Majer, 
S. Kumar, S. M. Girvin, and R. J. Schoelkopf. Nature, 431:162, 2004. 

[22] Brian David Josephson. Phys. Lett, 1:251, 1962. 

[23] Brian David Josephson. Advances in Physics, 14:419, 1964. 

[24] K.K. Likharev. Dynamics of Josephson Junctions and Curcuits. Gordon 
and Breach Science Publishers, New York, 1984. 

[25] Milton Abramowitz and Irene A. Stregun, editors. Handbook of Mathe- 
matical Functions. Dover Publishing Inc., New York, 1970. 



108 



[26] A.J. Leggett. Prog. Theo. Phys. Suppl, 69:80, 1980. 

[27] S.T. Ruggiero and D.A. Rudman, editors. Superconducting Devices. Aca- 
demic Press Inc., 1990. 

[28] Michael Thinkam. Introduction to Superconductivity. McGraw-Hill, Inc., 
NY, 1996. 

[29] R. Rouse, S. Han, and J.E. Lukens. Phys. Rev. Lett, 75:1614, 1995. 

[30] J.R. Friedman, V. Patel, W. Chen, S.K. Tolpygo, and J.E. Lukens. Na- 
ture, 406:42, 2000. 

[31] J.E. Mooij andd T.P. Orlando, L. Levitov, C.H. van der Wal L. Tian, 
and Seth Lloyd. Science, 285:1036, 1999. 

[32] I. Chiorescu, Y. Nakamura, C.J. P.M. Harmans, and J.E. Mooij. Science, 
299:1869, 2003. 

[33] R. Rajaraman. Solitons and instantons : an introduction to solitons and 
instantons in quantum field theory. North-Holland Pub. Co., 1982. 

[34] S. Nakajima. Progr.Theor.Phys., 20:948, 1958. 

[35] J.R. Senitzky. Phys. Rev., 119:670, 1960. 

[36] J. Prigogine and P. Resibois. Physica, 27:629, 1961. 

[37] R.P. Feynman and F.L. Vernon. Ann. of Phys., 24:118, 1963. 

[38] J. Zwanzig. J.Stat.Phys., 9:215, 1973. 

[39] O.A. Caldeira and A.J. Leggett. Phys. Rev. Lett, 46:211, 1981. 

[40] Ulrich Weiss. Quantum Dissipative Sysyems. World Scientific, 2 edition, 
1999. 

[41] W.K. Wootters. Phys. Rev. D, 23:357, 1981. 



109 



[42] D.V. Averin. Exploring the Quantun- Classical Frintier: Recent Advances 
in Macroscopics and Mesoscopic Quantum Phenomena. Nova Science, 
Hungtington, NY, 2003. also in cond-mat 0004364. 

[43] A.N. Korotkov and D.V. Averin. Phys.Rev B, 63:115403, 2001. 

[44] D.V. Averin. Phys.Rev. Lett., 88:207901, 2002. 

[45] A.N. Jordan and M. Biittiker. Phys.Rev. B, 71:125333, 2005. 

[46] B.P. Roe. Probability and Statistics in Experimental Physics. Springer, 
New York, 2 edition, 2001. 

[47] John W. Negle and Henri Orland. Quantum Many-Particle Systems. 
Westview Press, 1988. 

[48] K. Blum. Density matrix theory and applications. Plenum, NY, 1981. 

[49] M. Governale, M. Grifoni, and G. Schdn. Chem.Phys., 268:273, 2001. 

[50] M. Thorwart and P. Hanggi. Phys.Rev. A, 65:012309, 2002. 

[51] M.J. Storcz and F.K. Wilhelm. Phys.Rev. A, 67:042319, 2003. 

[52] J.R. Friedman and D.V. Averin. Phys.Rev. Lett, 88:050403, 2002. 

[53] Y. Nakamura, Yu.A. Pashkin, T. Yamamoto, and J.S. Tsai. Phys. Rev. 
Lett, 88:047901, 2002. 

[54] C. H. van de Wal, C. J. ter Haar, F. K. Wilhelm, R. N. Schoitern, 
C. Harmans, T. P. Orlando, S. Loyd, and J. E. Mooji. Science, 290:773, 
2000. 

[55] E. Paladino, L. Faoro, G. Falci, and R. Fazio. Phys. Rev. Lett, 88:228304, 
2002. 

[56] E. Paladino, L. Faoro, and G. Falci. Adv. Sol. State Phys., 43:747, 2003. 
[57] T. Itakura and Y. Tokura. Phys.Rev. B, 67:195320, 2003. 

110 



[58] Y.M. Galperin, B.L. Altshuler, and D.V. Shantsev. Fundamental Prob- 
lems of Mesoscopic Physics, pages 141-165. Kluwer Academic Publish- 
ers, The Netherlands, 2004. (also cond-mat/03 12490). 

[59] N.V. Prokof'ev and P.C.E. Stamp. Rep. Prog. Phys., 63:669, 2000. 

[60] J.L. Doob. Trans. Amer. Math. Soc, 52:37, 1942. 

[61] H.-P. Breuer and F. Petruccione. The theory of open quantum systems. 
Oxford Univ. Press, 2002. 

[62] Richard P. Feynman and A. R. Hibbs. Quantum Mechanics and Path 
Integrals. McGraw-Hill Companies, 1965. 

[63] K. Rabenstien, V.A. Sverdlov, and D.V. Averin. JRTP Lett, 79:646, 
2004. 

[64] K. Rabenstien, V.A. Sverdlov, and D.V. Averin. In H. Takayanagi and 
J. Nitta, editors, Realzing Controlable Quantum States - Mesoscopic Su- 
perconductivity and Spintronics in the Light of Quantum Computation. 
Worl Scientific, 2005. (also cond-mat/0401519). 

[65] Yu. Makhlin and A. Shnirman. JETP Lett., 78:497, 2003. (also cond- 
mat/0308297). 

[66] W.H. Press, S.A. Teukolsky, W.T. Vetterling, and B.T. Flannery. Nu- 
merical Recipes in C. Cambridge University Press, 1992. 

[67] R.Y. Rubenstin. Simulation and the Monte Carlo Method. John Wiley 
& Sons, 1981. 

[68] K.Y.R. Billah and M. Shinozuka. Phys. Rev. A, 42:7492, 1990. 

[69] E.O. Brigham. The Fast Fourier Transform. Prentice-Hall, 1974. 

[70] D.V. Averin, J.R. Friedman, and J.E. Lukens. Phys. Rev. B, 62:11802, 
2000. 



Ill 



[71] J. Dalibard, Y. Castin, and K. M0lmer. Phys. Rev. Lett, 68:580, 1992. 

[72] R. Dum, P. Zoller, and H. Putsch. Phys. Rev A, 45:4879, 1992. 

[73] C. W. Gardiner, A. S. Parkins, and P. Zoller. Phys. Rev. A, 46:4363, 
1992. 

[74] R. Dum, S. Parkins, P. Zoller, and H. Ritsch. Phys. Rev A, 46:4382, 1992. 

[75] H. Carmichael. An Open Systems Approach to Quantum Optics, Lecture 
Notes in Physics ml8. (Springer, Berlin, 1993. 

[76] J. A. Wheeler and W.H. Zurek, editors. Quantum Theory and Measure- 
ment. Princeton University Press, 1983. 

[77] P. Meystre and M.O. Scully, editors. Quantum Optics, Exsperimental 
Gravity, and Measurement Theory. Plenum, NY, 1983. 

[78] A.I. Miller, editor. Sixty Two Years of Uncertainty: Hystorical, Phy- 
losophical and Physical Inquires into the foundations of Quantum Me- 
chanics. Plenum, NY, 1990. 

[79] V.B. Braginsky and F.Ya. Khalili. Quantum Measurement. Cambridge, 
1992. 

[80] M.B. Mensky. Quantum Measurement and Decoherence: Models and Phe- 
nomenology. Kluwer, Dordtrecht, 2000. 

[81] I. Siddiqi, R. Vijay, F. Pierre, C.M.Wilson, M. Metcalfe, C. Rigetti, 
L. Frunzio, and M. H. Devoret. Rf-driven josephson bifurcation amplifier 
for quantum measurement. Phys. Rev. Lett., 93:207002, 2004. 

[82] W. Mao, D.V. Averin, R. Ruskov, and A.N. Korotkov. Phys. Rev. Lett, 
93:056803, 2004. (also cond-mat/0401484). 

[83] M. Field, C.G. Smith, M. Pepper, D.A. Ritchie, J.E.F. Frost, G.A.C. 
Jones, and D.G. Hasko. Phys. Rev. Lett., 70:1311, 1993. 



112 



[84] S.A. Gurvitz. Phys.Rev.B, 56, 1997. 

[85] M. Kataoka, J.B. Ford, G. Faini, D. Mailly, M.Y. Simmons, D.R. Mace, 
C.T. Liang, and D.A. Ritchie. Phys. Rev. Lett, 83:160, 1999. 

[86] A.N. Korotkov. Phys.Rev.B, 60:5737, 1999. 

[87] D. Sprinzak, M. Buks, E. Heiblum, and H. Shtrikamn. Phys. Rev. Lett., 
84:5820, 2000. 

[88] A.N. Korotkov. Quantum Noise in Mesoscopic Physics, page 205. Kul- 
wer, 2003. 

[89] CM. Caves, C.S. Thorne, R.W.P. Drever, V.D. Sandberg, and M Zim- 
mermann. Rev.Mod.Phys., 52:341, 1980. 

[90] V.B. Braginsky and F.Ya. Khalili. Rev.Mod.Phys., 68:1, 1996. 

[91] P. Grangier, J. A. Levenson, and J. P. Poizat. Nature (London), 396:537, 
1998. 

[92] S.Peil and G.Gabrielse. Phys. Rev. Lett, 83:1287, 1999. 
[93] W. Feller. Bull.Amer.Math.Soc., 51:800, 1945. 

[94] E. Il'lichev, N. Oukhanski, A. Izmalkov, Th. Wagner, M. Grajcar, H.G. 
Meyer, A.Yu. Smirnov, A. M. van den Brink, M.H.S. Amin, and A.M. 
Zagoskin. Phys. Rev. Lett, 91:097906, 2003. 

[95] A. Wallraff, D. I. Schuster, A. Blais, J. Mayer, M.H. Devoret, S.M. 
Girvin, and J.R. Schoelkopf. Phys. Rev. Lett, 95:060501, 2005. 

[96] A. Schmid. Phys. Rev. Lett, 51:1506, 1983. 

[97] Yu. Koval, A. Wallraff, M.V. Fistul, N. Thyssen H. Kohlstedt, and A.V. 
Ustinov. IEEE Trans. Appl.Supercond., 9:3957, 1999. 

[98] A. Shnirman, E. Ben- Jacob, and B. Malomed. Phys.Rev.B, 1997. 



113 



[99] M.B. Fogel, A.R. Bishop S.R. Trullinger, and J.A.Krumhansl. 
Phys. Rev. Lett, 36:1411, 1976. 

[100] D.W. McLaughlin and A.C. Scott. Phys. Rev. A, 18:1652, 1978. 

[101] A. Wallraff, A. Lukashenko, J. Lisefeld, A. Kemp, M.V. Fistul, Y. Koval, 
and A.V. Ustinov. Nature, 425:155, 2003. 

[102] D.V. Averin and E.V. Sukhorukov. Phys. Rev. Lett, 95:126803, 2005. 

[103] K.K. Likharev and V.K. Semenov. IEEE Trans. Appl.Supercond., 1:3, 
1991. 

[104] V.K. Semenov and D.V. Averin. IEEE Trans. Appl.Supercond., 13:960, 
2003. 

[105] A.M. Savin, J.P Pekola, D.V. Averin, and V.K. Semenov. cond- 
mat/0509318, 2005. 

[106] I.L. Aleiner, N.S. Wingreen, and Y. Meir. Phys. Rev. Lett., 79:3740, 1997. 

[107] Y. Levinson. Europhys.Lett, 39:299, 1997. 

[108] S. Pilgram and M. Biittiker. Phys. Rev. Lett., 89:200401, 2002. 

[109] A. A. Clerg, S.M Girvin, and A.D. Stone. Phys. Rev. B, 67:165324, 2003. 



114 



